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ABSTRACT 


^he  reliability  of  the  structure  can  be  calculated  by  combining  the 
non-uniform  stress  distribution  of  a  structure  from  finite  element  analysis 
with  the  statistical  strength  theory  under  two  dimensional  non- 
homogeneous,  uniaxial  stress  limitations.  The  specimen  models  with  and 
without  notch  are  the  sample  structures  to  illustrate  the  calculations.  The 
statistical  strength  of  the  structure  is  cast  in  the  standard  Weibull  form 

V  fej?  V 


characterized  by  the  structural  scale  parameter,  and  shape  parameter, 

a  U*v  o. 

5 cfe  which  are  functions  of  the  material  scale  and  shape  parameters  (w  and 

if*  ^ 

4.)  Therefore,  this  thesis  demonstrated  that  the  scale  parameter  of  the 


structure  remains  constant  for  different  load  magnitude,  only  the 


structural  geometry  and  loading  condition  cause  the  change  of  material 
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scale  parameter,  j  to  the  structural  scale  parameter,  tfE.  The  results  also 


shows  that  the  degree  of  uniformity  of  stress  within  structural  element 
and  the  linearity  of  the  load  magnitude  and  stress  distribution  within  the 
structure  affect  the  accuracy  for  the  calculation. 
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I.  INTRODUCTION 

A.  RELIABILITY  REQUIREMENT  FOR  COMPOSITE  STRUCTURES 

Reliability  theory  became  an  independent  scientific  discipline  at  the 
beginning  of  the  1950’s  under  the  influence  of  the  rapid  development  of 
radio  electronics,  computer,  devices,  and  rocketry.  Modern  radio  electronic 
apparatus  and  digital  mechanics  consist  of  a  very  large  quantity  of 
components.  If  the  failure  of  one  component  leads  to  the  failure  of  the 
unit  as  whole,  then  evidently  the  possibility  of  high  reliability  (freedom 
from  failures)  will  diminish  rapidly  as  the  unit  becomes  more  complex.  In 
this  connection,  problems  arise  on  the  prediction  of  the  reliability  or 
planned  units,  on  the  development  of  measures  to  increase  the  reliability, 
on  a  fundamental  for  reliability  testing  methods,  etc.  All  these  questions 
are  the  subject  of  reliability  theory. 

The  mechanical  strength  and  stiffness  of  a  structure  are  one  of  the 
aspects  of  reliability.  The  engineering  design  concept  is  not  only  solving 
the  stresses  and  strains  which  originate  in  structures  under  various 
external  imposed  conditions  but  also  trying  to  make  them  to  be  operated 
sufficiently  reliably  throughout  their  established  utilization  time. 
Therefore,  in  the  concluding  stages  of  an  engineering  design  structural 
mechanics  inevitably  comes  into  resolution  with  reliability  theory.  For 
critical  application  of  composite  as  an  aircraft  structure,  the  study  of 
reliability  in  composite  structures  is  required. 


B.  HISTORICAL  BACKGROUND 


Bullock  I  Ref.  I  ]  used  Weibull  theory  to  predict  the  strength  of  the 
composite  materials  between  different  laboratory  sample  configurations. 
He  investigated  the  strength  ratios  of  composite  materials  in  flexure  and  in 
tension  He  assumed  the  shape  parameter  (the  exponent  characterize 
flaw-density  that  determines  the  scatter  of  strength  of  the  material)  was 
not  changed  when  the  geometry  or  size  of  the  material  was  changed.  He 
used  the  same  shape  parameter  for  the  tensile  test  and  flexure  test  which 
based  on  homogeneous  state  of  stress  and  heterogeneous  state  of  stress 
respectively.  He  found,  however,  based  on  the  experimental  results,  the 
hypothesis  that  the  shape  parameter  was  constant  was  not  realistic. 
Changes  in  shape  parameter  were  in  fact  observed.  Whitney  and  Knight,  n 
[  Ref.  2  ]  used  a  statistical  strength  theory  based  on  Weibull  distribution  to 
explain  the  difference  between  unidirectional  tensile  data  generated  from  a 
flexure  test  and  a  standard  tensile  coupon  The  result  was  shown  a 
significantly  larger  variation  in  tensile  strength  versus  flexural  strength. 

Rosen  I  Ref.  3  ]  presented  a  theoretical  and  experimental  treatment  of 
the  failure  of  composites,  consisting  of  a  matrix  stiffened  by  uniaxially 
oriented  fibers,  when  subjected  to  a  uniaxial  tensile  load  parallel  to  the 
fiber  direction.  He  observed  that  a  portion  of  the  fiber  at  each  end  was  not 
fully  effective  in  resisting  the  applied  load  because  of  the  axial  load  was 
transmitted  by  shear  through  the  matrix  to  adjacent  fibers.  Basically,  the 
fibers  failed  as  a  result  of  a  statistically  distributed  flaws  or 
imperfections  and  composites  failed  as  a  result  of  a  statistical 
accumulation  of  such  flaws  over  a  given  region.  Therefore,  the  ineffective 


length  of  the  fiber  was  Introduced  in  order  to  consider  the  composite  to  be 
composed  of  a  series  of  layers  of  this  length.  Then  the  segment  of  a  fiber 
within  a  layer  might  be  considered  as  a  link  in  the  chain  that  contribute  the 
fiber.  Each  layer  was  then  a  bundle  of  such  links,  and  the  composite  as 
whole  can  be  modeled  by  series  of  such  bundles.  Weibult  theory  was  used 
to  define  a  statistical  distribution  of  flaw  or  imperfections  that  result  in 
fiber  failure  under  applied  stress.  The  statistical  accumulation  of  such 
flaws  within  a  composite  material  was  demonstrated  to  be  cause  of 
composites  failure.  He  concluded  that  the  stress  concentrations  occurred 
in  fiber  adjacent  to  the  fibers  break  were  considered  only  to  balance  the 
effect  of  the  variation  of  ineffective  length  which  affect  to  the  composite 
strength.  Knight,  C.E.  [  Ref.  4]  used  a  finite  element  stress  analysis  and 
statistical  strength  theory  to  assess  the  influence  of  the  stress 
concentration  on  the  ultimate  strength  calculation  for  composites  from 
the  Split-D  Test.  He  observed  that  the  geometry  changed  causes  the 
nonlinearity  and  the  boundary  condition  would  be  changed. 

The  Probabilistic  Statistic  Failure  of  Composite  Materials  [  Ref.  5  ] 
demonstrated  how  to  develop  a  theory  that  combined  the  statistics  of 
composite  material  failure  with  the  orthotropic  nature  of  composite 
materials.  The  effects  of  loading  history  and  the  probabilistic  location  of 
the  failure  for  a  composite  can  be  accounted  using  this  theory.  In  order  to 
improve  the  efficiency  of  composite  designs  for  structural  analysis,  a 
theory  that  combining  the  statistics  of  the  experimental  failure  data  with 
an  orthotropic  material  failure  response  was  developed.  The  analytical 
theory  developed  extent  deterministic  failure  laws  to  include  the 


statistical  scatter  observed  experimentally.  Tsai  and  Wu  [  Ref.  6  1  are  one 
of  the  investigators  who  described  a  general  theory  for  the  failure  of 
composite  materials  as  the  strength  tensor  theory  but  determining  is 
complicated  by  wide  range  in  failure  properties  of  composite  materials  and 
scatter  in  experimental  results.  The  Weibuli  weakest  link  theory  was 
developed  and  described  a  description  of  the  spatial  distribution  of  the 
failure  locations  with  a  sample  calculation  performed  for  the  case  of  a 
uniformly  loaded  tensile  specimen.  Three  point  bending  specimens  were 
used  to  predict  failure  location  by  applying  the  basic  weakest  link  theory. 
Four  point  bending  specimens  were  used  to  illustrate  the  effect  of  loading 
history  which  were  evaluated  by  applying  numerical  integration  to  the 
probability  of  failure  density  distribution  curves.  The  compressive  failure 
load  was  demonstrated  that  it  is  not  deterministic  failure  and  also  having 
some  scatter  which  is  independent  of  tensile  scatter  so  that  the  Weibuli 
failure  distribution  for  tension  and  compression  are  different. 

In  ordinary  structures  having  many  stress  components  active,  a 
decision  on  the  probability  for  survival  must  be  based  on  statistical 
descriptions  that  admit  multiaxial  stress.  The  development  of  combining 
weakest  link  theory  with  the  strength  tensor  representation  is  required. 

C.  OBJECTIVE  AND  SCOPE 

The  driving  motivation  for  this  research  is  that  many  critical 
components  of  modern  aircraft  structure  are  made  of  composite  materials. 
However,  the  reliability  design  is  necessary  for  engineering  designs 
requires  not  only  the  stress  analysis  but  also  statistical  theory  of  failure 


strength  of  composites,  The  latter  motivated  the  investigation  of 
statistical  theory  for  composites  failure  strength  where  are  summarized  in 
the  last  sectioa  Any  aircraft  structural  components  inevitability  have  the 
non-uniform  stress  in  the  neighborhood  of  edge-notches  and  holes.  Size 
effect  is  the  one  that  is  importance  for  the  composite  study  because  of  the 
difficulty  and  the  economy  that  one  could  not  experiment  the  several  large 
samples  to  predict  the  reliability  of  the  components.  Even  the 
experimental  test  samples  can  be  made,  the  difficulty  to  identify  the 
failure  location  is  also  impossible.  Finite  element  method  and  the 
statistical  theory  can  be  combined  to  explore  the  possibility  of  prediction 
both  the  reliability  of  the  structure  and  the  failure  location.  Basically, 
finite  element  method  can  be  used  to  calculate  the  stress  distribution  that 
occurs  in  the  structure  under  an  arbitrary  loading  boundary  condition.  The 
structure  is  divided  into  elements  be  sufficiently  small  that  the  stress  is 
uniform.  This  solution  of  physically  continuous  system  can  then  be 
replaced  by  the  solution  of  the  discrete  system  of  the  elements.  The 
statistic  Weibull  theory  can  be  used  to  model  composite  materials 
composing  of  a  series  of  bundles  which  are  in  turn  composited  series  of 
fiber  links.  Another  assumption  for  the  statistical  Weibull  theory  is  the 
local  stress  in  each  link  must  be  uniform.  For  the  connection  to  the  two 
dimensional  finite  element  method,  instead  of  using  link,  the  element  will 
be  introduced.  For  one  dimensional  case,  the  link  can  be  used  but  in  two 
dimensional  case,  the  element  has  a  width.  The  operation  of  size  effect 
formulation  based  on  Weibull  theory  and  finite  element  method  will  be 
presented  in  next  chapter. 


The  scope  of  this  research  limits  to  the  treatment  of  statistical 
strength  of  composite  structures  which  are  two  dimensional,  subjected  to 
nonhomogeneous  uniaxial  stress  including  size  effect.  The  specimen 
models  are  plate  with  and  without  notch  subjected  to  the  tensile 
displacement  loading  case.  The  finite  element  code  name  NIKE2D  [  Ref.  7  ] 
is  used.  NIKE2D  is  a  two  dimensional  finite  element  program  which  is 
developed  by  Hallquist,  J.O.  [  Ref.  7  L  It  can  be  used  to  calculate  linear  and 
non-linear  finite  element  problems,  in  this  thesis,  only  the  linear  portion 
is  utilized.  MAZE  [  REf.  8  ]  computer  program  is  used  to  generate  specimen 
model,  grid,  define  boundary  conditions,  and  also  create  the  input  file  that 
can  be  immediately  use  for  NIKE2D  program.  ORION  [Ref.  9]  is  the  post 
processing  of  NIKE2D  which  can  be  used  to  create  contour  curves  of  stress, 
displacement,  etc.  The  observation  of  uniform  stress  in  each  element  will 
be  made.  The  observation  of  using  half  of  the  entire  model  by  symmetrical 
will  be  made,  which  bases  on  the  idea  that  the  model  can  be  divided  into 
smaller  elements  when  compare  to  the  entire  model.  The  linear 
relationship  between  external  loads  and  the  local  stress  in  each  element 
will  be  verified.  The  plate  without  notch  is  used  as  bench  mark  to  compare 
the  reliability  results  to  the  plate  with  notch.  The  specimen  model 
dimension  bases  on  specimen  dimension  which  using  in  the  shear  test 
experiment.  The  shear  model  that  involves  with  the  muitiaxiai  stress 
which  can  not  yet  be  treated  in  this  thesis.  Future  work  requires 
generalization  of  current  results  to  muitiaxiai  stress.  When  such  extension 
is  made,  realistic  prediction  of  composite  structures  will  then  be  possible. 


A.  REVIEW  OF  WE1BULL  WEAKEST  LINK  FAILURE  THEORY 

consider  an  infinitesimal  volume  element  SV  at  point  x  and  subjected 

to  stress  state  Ojj.  The  probability  of  failure  Sf  of  volume  element  SV 


is  assumed  to  be  given  by: 


Sf  =  *  [  Ojj  [x]  ]  SV 


(  2.1  ) 


where; 


¥  =  failure  parameter  which  is  a  function  of  Ojj  which  is  a 


function  of  position,  x. 

It  is  assumed  there  are  only  two  possible  states  for  the  volume 
element,  SV  ,  either  the  volume  element  survives  the  loads  or  fails.  From 
this  assumption  the  probability  for  survival  or  reliability  is  simply: 

SR  =  1  -  Sf  =  1  -  "H  Ojj  lx]  ]  SV  (2.2) 


The  basic  assumption  in  weakest  link  theory  is:  If  any  volume,  SV  ,  in 
the  total  structure  volume,  V  ,  fails,  then  the  entire  structure  occupying 
volume  V  fails. 

Consider  the  structural  volume,  V,  to  be  divided  into  N  volume 
elements  SVk  ,  with  each  of  the  volume  elements  located  at  point  xk  . 


Then  the  total  volume  is  given  by: 


V  =  liml  SVk  =  Jv  dV;  k  =  1,  2...  N 


(2.3) 


The  probability  for  the  volume  V  to  survive,  R  ,  under  weakest  link 
theory,  is  given  by: 

R  =  lim  IT  SRk  =  lim  TT  { 1  -  ♦[  Ojj[xk1  ]  8Vk }  (  2.4  ) 

where:  k  =  l,  2  . . .  N ,  and 

N-*oo 

Equation  2.4  can  be  equivalently  written  as: 

In  R  =  lim  E  In  {  I  -  Ojj[xkl  SVk  }  (2.5) 

where:  k  =  I,  2  . . .  N ,  and 

N-  00 

The  quantity  *.8Vk  is  small  and  from  the  Taylor  series  expansion 

In  (1  ♦  x)  s  x  for  x  «  1 

Equation  2.5  becomes: 

In  R  =  lim  E  -{ aV]  [xkl  J  SVk }  (  2.6  ) 

or  R  =  exp  -{  Jv  ¥[  cfy  [xkl  ]  dV  }  (  2.7  ) 

where:  k  =  I,  2  . . .  N  ,  and 

N  -*<*> 

Equation  2.7  is  the  basic  equation  derived  under  the  assumption  of 

weakest  link  theory.  Note  that  the  stress  distribution  need  not  to  be 
uniform;  the  material  can  be  generally  anisotropic  and  the  material  can  be 
non-homogeneous.  The  equation  is  limited  in  that  it  describes  only  the 
reliability  of  the  entire  structure  and  provides  only  a  quantitative 
description  of  where  the  structure  will  fail.  The  reliability  becomes: 

R  =  exp  -  { Oii/B  )* 


where*. 


f  m  S 


•  L  *  -'V  r*  .  w  .  *  «  .  '  I 


<x  =  shape  parameter,  and 
3  =  scale  parameter 
In  the  next  section,  basic  theory  will  be  used  to  obtain  the 
formulation  of  size-effect  for  one  dimensional,  non-homogeneous, 
uniaxial  stress. 

B.  SIZE-EFFECT  FOR  NON-HOMOGENEOUS  STRESS,  ONE  DIMENSIONAL, 

UNIAXIAL  STRESS 

The  value  of  <xLj  arxl  #LI  are  estimatBCl  by  the  experimental 
procedure  from  the  given  n  numbers  of  sample  length  LI .  For  a  sample 
length  L  ,  the  evaluation  of  <xL  and  3|_  can  be  implemented  according  to 
the  following  procedure. 

There  are  j  links  in  length  L  (Figure  2.1-a)  each  links  are  not 
necessarily  of  equal  length,  but  the  stress  in  each  link  is  homogeneous, ie 
the  stress  oxp  within  the  p*h  link  satisfies  the  following  condition: 

xoxp  -  {/|p  xoxdx}/lp  $  e  (  2.8  ) 

where:  €  can  be  any  arbitrarily  small  value 

Ip  =  length  of  pth  link 

P  =  I.  2  .  . .  j 

The  reliability  for  each  link  is: 

Rp(tfxp)  =  I  -  f(cxp)  =  exp  -(axp/£p)°<p  (  2.9  ) 
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w./.v.v: 


.  v-  • 


There  are  (lp/1)  unit  links  within  lp  ,  the  reliability  of  each  unit  link 
within  lp  is: 

Rp(oxP)  =  exp  -(<JxP/$  )“  (  2.10  ) 

where:  «  and  3  for  unit  link  are  material  constants 

The  reliability  for  each  link  lp  is: 

Rp(oxP)  =  TT  {exp  -(axP/$)“}m  ;  m  =  I,  2  .  .  .  (Ip/1) 

=  exp  {-(lp/l)(a 

where:  (lp/l)  is  the  number  of  unit  link  in  lp 

The  reliability  for  then  entire  length  L  is: 

R|_ =  Rj  RsC^)  04(^4)  •  •  •  Rj(Oj) 

=  exp  {-(li/1)(d  1/3)“}  exp  {-(I2/IX0  2/3)“)-  •  •  exp  {-(Ij/lXo  1/3)“] 

=  exp  -{(li/0(cf  i/3)<x  ♦  O2/1X0  2/3)<x  4  •  •  •  ♦  0j/0(cf  j/3)<x) 

Rl  =  exp  -{ Z  Klp/lXcr  xp/3)<xl )  ;  p  =  1.  2  . . .  j  (  2.11  ) 

Under  uniform  stress  condition  within  an  element,  let  op  is  the 

element  stress  which  is  constant  within  lpth  length  element  and  op  is 

calculable  from  a  given  boundary  condition  (B.C.)  characterizable  by  {P}, 
(eg,  a  load  vector).  The  stress  due  to  the  same  B.C.  with  a  different 
magnitude  P  is: 

oplP]  =  dp[{P)]  (P/CP))  (  2.12  ) 

where:  P  can  be  of  any  arbitrary  magnitude  in  linear  system,  or 

P  is  of  neighboring  magnitude  of  local  linear  system. 


opl{P)l  is  expressible  from  structural  mechanics  calculations, 
(eg,  a  finite  element  method.) 

In  order  to  reduce  the  value  of  oxp  which  is  a  function  of  P  in 
equation  2.11  to  a  single  random  variable  {P}  ,  the  substitution  of  oxp[P] 
by  0XPUP}1  is  made  within  the  same  manner  as  the  expression  in  equation 
2.12.  Thus,  the  equation  2.11  becomes: 

Rl  =  exp  -{ I  C(lp/lXcr  xpl{PMP})*(P/3)“] )  ;  p  =  I,  2  . . .  j 

Expressing  the  above  equation  into  the  Standard  Weibull  form: 

Rl  =  exp  -{ (P/3L)°<L  } 

Then:  -  cx 

flL  =  fl/  ( z  [(Ip/lXo  K^fPll/lP})*!  }IA*  ;  p  =  I,  2  . . .  j 

Let:  F  =  ( z  [(lp/l)  i/k  (axp[(P)l/(P))l“  )-|/K  :  p  =  I,  2  . . .  j 

Thus:  3j_  =  3F 

That  is,  given  a  structure  of  length  L  subjected  to  boundary  condition 
parameter  (P)  (  which  may  be  of  a  single  or  multiple  load  Pi  ,  P2,  .  .  . 
vectors  characterizable  by  a  single  parameter  {P}  )  which  give  rise  to  a 
unidirectional  stress  op  computable  from  the  methods  of  mechanics  (  eg, 

finite  element  analysis  ).  Ip  are  element  lengths  which  are  segmented  as 
small  as  necessary  to  assume  the  stress  to  be  uniform.  Upon  sampling  n 
structures  measuring  Pn  failure  loads. 


«L ,  fiL  are  shape  and  scale  parameters  ( in  dimension  or  load 
estimated  from  D(Pn)  and  structural  size  L  . 

<x  ,  3  are  shape  and  scale  parameters  ( in  dimension  of  stress  which 
are  geometric  independent  material  parameters  ). 

F  is  a  conversion  factor  calculable  for  structural  (  finite  element  ) 
analysis. 

For  special  case,  homogeneous,  one  dimensional  uniaxial  stress: 

From:  F  =  { £  [(lp/l) ]/<*  (a/KPMP})]*  )'Uo(  ;  p  =  1,  2  . . .  j 

For  axp  =  a  and  {P}  =  Pn: 

F  =  [  (a/Pn)<x  (  S  (lp/l)  l/“  )“  rl/<x  ;  p  =  I,  2  ...  j 

=  I  (<J/Pn)  {  Z  (lp/1)  ,/k  )  l'1  ;  p  =  I.  2  . .  .  j 

=  (<J/Pn)-'  (L/1)  -|/a 

Thus: 

Rl  =  exp  -{ [  (P/BJAa/Pn)’1  /(L/1)  -,/k  /" ) 

=  exp  -{ (L/1)  [  (a/flXP/Pn)  ]« ) 

This  recovers  the  well  known  classical  result  demonstrates  the 
consistency  of  the  current  formulation  is  true. 

For  each  given  non-homogeneous  boundary  condition,  a  corresponding 
correction  factor  F  can  be  calculated  by  finite  element  such  that  the  size 
effect  can  be  predicted.  By  weakest  link  in  length  X  ;  it  is  inferred  that 
the  plane  of  facture  is  perpendicular  to  X  .  That  means  the  failure  plane 


(4  |L  iL  »i  tta  aL  ila'il  %)  «<i  • 


has  no  thickness.  Since  the  one  dimension  specimen  has  no  width,  fracture 
width  is  not  a  parameter. 

C.  SIZE-EFFECT  FOR  NON-HOMOGENEOUS  STRESS,  TWO  DIMENSIONAL, 

UNIAXIAL  STRESS 

This  is  the  extension  of  one  dimensional,  uniaxial  stress  to  two 
dimensional  uniaxial  stress.  The  definition  that  uniaxial  stress  implies 
failure  has  zero  dimension  in  X  direction  is  adopted  for  one  dimensional, 
uniaxial  stress  ,  which  has  been  described  in  the  previous  section.  For  two 
dimensional,  uniaxial  stress,  failure  automatically  requires  definition  of 
the  failure  plane  because  the  specimen  has  dimension  in  Y  direction 
(Figure  2.1-b),  and  it  is  no  longer  conceptually  a  point  as  in  the  one 
dimensional  case.  In  the  similarly  procedure  as  in  one  dimensional  case, 
the  specimen  is  divided  into  sufficiently  small  areas  within  which  the 
stress  in  direction  of  force  is  assumed  to  be  uniform.  It  is  seen  that  not 
only  the  length  of  the  specimen  is  divided  into  small  element  lengths  but 
also  the  width  must  be  divided  into  small  element  widths  in  order  to 
provide  element  areas  which  are  satisfies  uniform  stress  assumption.  Thus: 


xoxP  -  jpJ  xoxdx/lp  $  €x  for  all  lf 


where; 


P  =  I.  2  . . .  j 


y°xP  "  wqJ  yaxdy/wq  *  €y  for  311  wq 
where:  q  =1,  2  . . .  k 

Weakest  link  in  X  direction  is  defined  by  Eq.  2.8  and  unit  length  in  X 
direction  is  defined  by  Eq.  2.9. 


"• . \  *•  \ 
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Weakest  link  in  Y  direction  is: 


Rqloxpq1  =  exp  -{cJx^/^pqJ^q 

where:  oxp<l  is  stress  in  the  q*h  element,  and  q  =  I,  2  . . .  k 

This  means  that  if  any  one  of  the  q*h  fail,  all  elements  fail.  However, 
the  unit  (metric)  width  can  not  be  arbitrarily  defined.  It  is  visualized  that 
within  each  metric  area,  there  is  one  domineering  crack  (largest  crack 
perpendicular  to  Ox).  When  subjected  to  oxp<J,  the  domineering  crack 
extends  without  bound.  From  fracture  mechanics: 
kc  =  o/a^ 

where:  kc  is  a  material  constant  and  can  be  independently  measured. 

ac  is  weakest  link  in  Y  direction,  which  can  be  used  as  the  unit 
metric  width. 

Thus:  acp<J  =  (kj/o^)2  for  element  pqth 

As  well  as  the  division  of  length  element  lp  into  unit  length,  the 

division  of  width  element  is  necessary.  In  stead  of  unit  length  I  ,  acp<*  is 

used  in  order  to  indicate  the  critical  width  in  each  element,  thus,  the 
reliability  of  entire  element  is: 

Rpql°xP<^  =  '<  Opq/D(w  pq/acpqXa1(pq/3)0< ) 

=  exp  -(  (Ipq/lXw  M/(k<70„pq)2](axp<l/3)°! ) 

=  exp  -( (lpq/l)(w  p(]/kc2)(<Jx>><!)2*lx/3°< ) 
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where:  pq  =  1,  2  .  .  .  jk  (  subscript  is  changed  from  two  dimensional 

array  into  one  dimensional  array.) 

The  reliability  of  entire  specimen  (which  is  the  structure)  is: 

Re  =  R,  R2  R3  . . .  Rjk 

=  exp  -{  Z  (lpq/1)(w  pq/K2)(oxP^)2+0</30<  } ,  pq  =  1.  2  . . .  jk 
From  the  linear  relationship  between  external  load  and  local  stresses: 
oxP<![P]  =  a^[{P)]  (P/(P» 

where;  axp^{P}]  is  element  stress  which  is  calculated  under  certain 
loading  boundary  condition  while  the  magnitude  of  load  is  {P}. 

oxP<HP1  is  element  stress  which  subjected  to  the  same 
boundary  condition  as  above  but  different  in  magnitude,  the  magnitude  is  P. 

Thus:  Rj  =  sup  -(  E  (lpq/l)(w  pq/kc2X(P/fP})Vq[(P)])2*<x/J0<  } 
where:  pq  =  1,  2  . . .  jk 

Express  it  in  term  of  the  Standard  Weibull  for  the  structure,  thus: 

RE  =  exp  -{ (P/3E)°<E  (  2.13  ) 

where:  «E  =  «  ♦  2 

%  =  z  Opq/lXw  pq/k^Xo^P)!)-**2  )(-'A«*2)) 

where:  pq  =  I,  2  . . .  jk  (  2.14  ) 

The  values  of «  and  3  of  material  are  altered  by  the  boundary  condition 
and  the  geometry  of  the  specimen  (or  structure)  changed  but  they  do  not 
depend  on  the  changed  of  the  magnitude  of  the  load.  Therefore,  the  shape 
parameter,  <x  and  the  scale  parameter,  3  can  be  solved  for  any  given 


structural  geometry  and  boundary  condition  The  f mite  element  method  can 
be  used  to  calculate  the  stress  distribution  of  the  specimen  by  dividing  the 
specimen  into  small  elements  within  which  the  stress  element  is  uniform. 

For  the  simple  limiting  2-D  case,  the  specimen  without  notch,  the 
stress  distribution  for  entire  specimen  is  uniform.  Thus: 

cxp<^  =  ox  =  constant 

SE  =  «P}/<J„)  s  [lp<)Wp()|  j(-l/«+2) 

=  (P)  |;c(2/o<e)  a(1/«e) 


where:  a  =  Z  l  lpqWpq  1 

<xE  =  ot+2 


,  pq  =  1,  2  . . .  jk 


Thus:  Re  =  exp  -{ (A/kc2)  (P/{P})  (ox/3<x)  )*£ 


(  2.15  ) 


.*•  ”•  l.  "  ■<  .*•  k  k'p^A  .  '■  S  , 


A.  SPECIMEN  MODEL  CONFIGURATION  AND  COORDINATE  SYSTEM 

The  specimen  model  has  1.99  *  in  length  and  0.968*  in  width  is  align 
in  the  y-z  plane,  (cartesian  right  hand  coordinate  system  is  used).  The 
configuration  and  the  coordinate  of  the  specimen  without  notched  and  the 
plate  with  90  degree  notch  configuration  are  shown  in  Figure  3.1.  The 
fiber  direction  which  makes  the  orthotropic  material  having  large  elastic 
modulus  (  Ejj  )  is  considered  to  be  along  the  y-axis  (Figure  3.2). 

B.  FINITE  ELEMENT  PROCEDURE 
1.  Material  Constants 

The  finite  element  program  requires  defining  the  material 
constants  [  Ref.  7  ]  for  orthotropic  elastic  to  be  used  in  stiffness  matrix 
calculatioa  The  material  constants  of  typical  graphite  epoxy  composite 
which  used  in  this  thesis  are: 

Ejj  =  18.65  Msi  (Mega  Pound  per  square  inch) 

E22  =  1-21  Msi 
E,t  =  1.21  Msi 


Gj2  1  0.84  Msi 

2.  Boundaru  Conditions 

The  specimen  model  is  subjected  to  the  displacement  boundary 
condition  (B.C.)  in  order  to  be  able  to  use  the  entire  specimen  model  in  the 
calculation.  An  original  problem  the  specimen  model  is  subjected  to  the 
displacement  B.C.  value  ut  (Figure  3.3)  on  the  right  hand  side  in  the 
positive  y-direction  and  the  left  hand  side  is  subjected  to  the  displacement 
B.C.  value  u1  in  opposite  direction.  In  order  to  fix  the  specimen  model 

from  rigid  body  displacement,  the  left  side  of  the  model  must  be  constraint 
in  the  direction  of  load.  Therefore,  the  displacement  value  on  the  right 
side  of  the  model  becomes  two  time  of  Uj.  The  transform  of  the 

displacement  B.C.  is  shown  in  Figure  3.4. 

3.  Data  Planning  and  Purposes 

NIKE2D  is  two  dimensional  axis  symmetry  finite  element  code 
that  the  thickness  is  assumed  to  be  infinite.  The  plane  strain  analysis  is 
used  by  assuming  no  displacement  along  x-axis  and  no  rotation  around  y  and 
z  axis,  which  it  is  satisfy  for  this  purpose.  There  are  eight  finite  element 
program  data  output  to  be  discussed  in  this  thesis: 

3.1  Specimen  model  without  notch  subjected  to  the 
displacement  B.C.  value  0.01” 

3.2  Same  specimen  as  3.1  but  the  displacement  B.C.  value  is 
0.001".  The  purpose  of  these  two  cases  is  to  demonstrate  that  there  is  no 
stress  concentration  along  displacement  direction  and  the  linear 
relationship  between  displacement  and  local  stress. 
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3.3  Specimen  model  with  notch  subjected  to  the  displacement 

B. C.  value  0.01*  and  the  stress  calculation  bases  on  four-point  integration 
(Ref.  7 )  in  order  to  observe  the  uniform  stress  in  each  element. 

3.4  The  same  as  3.3  but  the  calculation  bases  on  one-point 
integration  at  the  centroid  of  each  element. 

3.5  Half  of  the  specimen  model  with  notch  subjected  to  the 
displacement  B.C.  value  0.01*  and  four-point  integration  is  used.  This  case 
will  be  used  to  compare  with  the  entire  specimen  model  in  order  to  insure 
that  the  half  specimen  can  be  used  because  of  the  symmetrical 
configuratioa 

3.6  Same  as  3.5,  but  one-point  integration  is  used. 

3.7  Same  as  3.6,  but  the  displacement  B.C.  value  is  0.004. 

3.8  Same  as  3.6,  but  the  displacement  B.C.  value  is  0.001. 
One-point  integration  results  will  be  used  to  be  input  data  for 

statistical  program. 

C.  STATISTICAL  PROCEDURE 

After  obtaining  the  results  from  the  finite  element  program, 
four -point  integration  results  have  to  be  considered  to  insure  that  the 
tensile  stress  in  each  element  is  uniform  or  almost  uniform.  If  they  are 
not  uniform,  the  size  of  the  element  need  to  be  reduced.  However,  due  to 
computer  memory  limitation,  absolute  uniform  stress  may  not  be 
attainable.  One  should  realize  that  what  spatial  domain  that  the  element 
size  should  be  reduced  and  what  area  the  element  need  not  be  changed.  The 
symmetrical  of  the  model  can  be  considered  to  be  used  only  half  of  the 
entire  model  in  order  to  increase  the  number  of  elements. 


The  statistical  size  effect  in  two  dimensional,  non-homogeneous, 
uniaxial  stress  program  in  fortran  77  is  written  to  be  used  to  calculate  the 
reliability  of  the  specimen  models  by  using  the  finite  element  results  and 
shape  parameter, «  and  scale  parameter,  3  from  several  numbers  of  certain 
sample  size.  There  is  one  program  including:  I)  the  specimen  model 
without  notch,  2)  the  entire  specimen  model  with  notch,  and  3)  the  half 
(top  side)  specimen  model  with  notch. 


A.  EXAMINATION  FOR  UNIFORMITY  OF  UNIFORM  STRESS  WITHIN 

ELEMENT 

The  four-point  integration  are  obtained  in  entire  and  half  specimen 
models  subjected  to  the  displacement  boundary  condition  (B.C.)  value  0.01. 
For  the  full  specimen  model,  the  element  number  544  (shade  in  Figure  4.1) 
is  considered,  the  tensile  stress  values  can  be  read  in  figure  4.2  .where  the 
stresses  of  this  element  are  tabulated  under  rows  2173,  2174,  2175,  and 
2176  (relatable  to  the  element  by  dividing  the  last  number  by  4  recovering 
544)  and  the  tensile  stresses  are  0.03663  Msi,  0.03738  Msi,  0.01933  Msi, 
and  0.01851  Msi  (Figure  4.2)  respectively,  which  demonstrated  that  the 
stress  within  element  is  not  exactly  uniform  but  the  differences  are  small. 
The  element  number  of  the  half  specimen  model  at  the  same  position  as 
element  number  544  of  the  full  specimen  is  442  (shade  in  Figure  4.3)  and 
the  tensile  stress  can  be  read  in  Figure  4.4,  where  the  row  numbers  are 
1765,  1766,  1767,  and  1768  and  the  tensile  stresses  are  0.03667  Msi, 
0.03742  Msi,  0.01937  Msi,  and  0.01861  Msi  respectively,  again  they  are  not 
exactly  uniform  but  the  differences  are  small.  The  element  number  192  in 
Figure  4.5  (full  specimen  model)  is  outside  the  notch  region,  the  tensile 
stress  values  can  be  read  in  Figure  4.6,  where  the  row  numbers  are  765, 
766,  767,  and  768  and  the  tensile  stresses  are  0.06891  Msi,  0.06888  Msi, 
0.06852  Msi,  and  0.06855  Msi  respectively  ,  which  shown  that  the  stress 
within  element  is  not  uniform  but  the  differences  are  more  uniform  than 
element  in  the  notch  region.  The  element  number  of  the  half  specimen 


model  at  the  same  position  as  element  number  192  of  the  full  specimen  is 
230  (Figure  4.7)  and  the  tensile  stress  can  be  read  in  Figure  4.8,  where 
the  row  numbers  are  917,  918,  919,  and  920  and  the  tensile  stresses  are 
0.06906  Msi,  0.06905  flsi,  0.06871  flsi,  and  0.06872  flsi,  while  they  are  not 
exactly  uniform  but  the  differences  are  also  smaller  than  the  element  in  the 
notch  regioa  in  principle  the  stress  in  each  element  may  approach 
uniformity  by  decreasing  the  size  of  the  element.  This  will  increase  the 
number  of  the  elements  and  is  limited  only  by  the  computer  memory  and 
computation  time.  The  observation  of  symmetrical  model  is  obvious  from 
the  stress  values  of  full  specimen  model  and  half  specimen  model. 

B.  OBSERVATION  OF  LINEAR  RELATIONSHIP  BETWEEN  EXTERNAL  LOAD 

AND  LOCAL  TENSILE  STRESS 

The  tensile  stress  of  the  element  number  460  (Figure  4.3)  of  half 
specimen  model  which  subjected  to  the  displacement  B.C.  value  0.001  is 
0.002475  Msi  (Figure  4.9)  and  the  tensile  stress  of  the  same  element 
which  subjected  to  the  displacement  B.C.  value  0.004  is  0.0112  Msi  (Figure 
4.10).  The  ratio  of  tensile  stress  is  almost  the  same  as  the  ratio  of 
displacement,  the  different  value  bases  on  the  non-uniform  of  the  stress 
around  the  notch  region  and  can  be  accepted.  When  the  element  number  230 
(Figure  4.7)  of  the  same  model,  which  is  far  from  the  notched  region  the 
ratio  of  the  stresses  (Figure  4.11,  displacement  B.C.  value  0.001  and 
Figure  4.12,  displacement  B.C.  value  0.004)  and  the  ratio  of  the 
displacement  are  almost  the  same  and  the  difference  is  small  compare  to 
the  different  value  for  the  element  in  the  notch  regioa  Therefore,  the 
linear  relationship  between  external  load  and  local  stress  can  be  accepted 


whether  the  element  is  in  notch  region  or  not.  Since  the  element  is  small 
enough  to  be  satisfy  the  uniform  stress  within  element,  then  the  uniformity 
and  the  linearity  can  be  attained. 

C.  STATISTICAL  SIZE  EFFECT  RESULTS  DISCUSSION 

Equation  2.13  and  2.14,  calculated  the  reliability  of  the  specimen 
model  based  on  statistical  size  effect  in  two  dimensional,  non- 
homogeneous,  uniaxial  stress.  For  application  in  these  equation,  the 
displacement  boundary  condition  need  to  be  converted  to  the  average  stress 
along  the  boundary.  The  definition  of  the  average  stress  in  two  dimension 
is; 

(P}=  { £  [wpqoxp<l] }/  { S  wpq } 

where:  pq  (element  number)  =  186,  187,  188,  189,  190,  200,  205,  210, 
215,  220,  225,  and  230  for  half  specimen  model  (Figure  4.7),  and 

pq  =  163,  164,  165,  166,  167,  168,  171,  174,  177,  180,  183,  186, 
189,  192,  465,  468,  471,  474,  477,  480,  483,  and  486  for  entire  specimen 
(Figure  4.5  and  Figure  4.13). 

The  values  of  stress  in  each  element  are  shown  in  Table  4.1  and  Table 
4.2  respectively. 

The  statistical  strength  of  a  structure  is  characterized  by  the  scale 
parameter  3f  in  equation  2.13  (  same  as  BET  A  A  in  computer  program. 

Appendix  A).  A  change  of  the  3e  value  is  therefore  related  to  the  strength 

change  of  the  structure.  The  values  of  3e  (equation  2.14)  for  plate  with 

notch,  where  the  values  of  <x,  3.  and  kc  are  10,  0.1  Msi,  and  0.001 


0.001  respectively  are  shown  in  Table  4.  3- A.  Table  4.  3-A  also  shows 
that  the  variation  of  the  values  of  0E  which  are  calculated  from  different 
values  of  external  load  (PBAR)  from  different  values  of  displacement 
boundary  condition  (u').  The  values  of  JJE  which  calculated  from  half 

specimen  model  seem  to  be  higher  than  %  value  from  the  full  specimen 

model.  The  reason  is  that  the  number  of  elements  for  full  specimen  is  less 
than  the  number  of  elements  for  half  specimen  model,  therefore,  the 
accuracy  for  full  specimen  model  could  be  lower  than  that  for  half 
specimen  model.  For  the  purpose  of  reliability  calculation  for  the 
structure,  one  can  average  the  0E  values  from  one  set  of  finite  element 
data  output  from  different  values  of  load  under  the  same  boundary 
condition,  therefore, for  instance  the  average  value  of  0E  is  approximately 

0.030.  Table  4. 3-B  shows  the  values  of  JJE,  where  the  values  of  <x,  and  kc 


are  changed  to  20,  and  0.0001  respectively.  The  0E  values  are  reduced  but 
the  difference  of  the  average  0E  values  is  not  large,  because  of  the 
compensation  of  higher  <x  and  lower  kc  values.  This  means  that  the 

changes  of «  and  kc  values  do  not  affect  much  to  the  3E  value  but  <x  affects 
the  reliability  of  the  structure  (equation  2.13).  Ail  the  reliability 
calculation  results  for  different  values  of  0E  .where  the  specimen  models 
are  operated  in  the  same  loading  range  are  shown  in  Appendix  B. 


In  the  case  of  plate  without  notch  that  shown  in  Table  4.  3-C  and 
Table  4.  3-D,  the  values  of  from  different  values  of  displacement 
boundary  condition  are  exactly  the  same.  The  reason  is  that  the  uniformity 
and  the  linearity  affect  the  accuracy  of  the  calculations.  The  value  of 

from  the  plate  without  notch  (0.0439)  is  higher  than  from  the  plate  with 
notch  (0.028)  under  the  same  boundary  condition  (u'  =  0.01),  which  is  a 
reasonable  result. 

In  addition,  the  results  show  that  the  material  scale  parameter,  3  is 
converted  to  3f  by  the  geometry  of  the  structure  and  the  loading  boundary 

condition,  but  it  does  not  depend  on  magnitude  of  the  load.  The  results  of 
the  specimen  model  without  notch  confirmed  that  the  uniform  stress 
assumption  is  valid  and  known  result  can  be  recovered.  Between  results  of 
notched  specimen  model  we  showed  that  3g  remains  constant  for  different 

load  magnitude  because  the  uniformity  and  the  linearity  conditions  are 
attained.  As  expected  we  confirmed  the  notched  model  is  less  reliability 
than  the  plain  model  under  the  same  loading  condition.  Given  available 
experimental  data  of  a  composite  material,  we  can  evaluate  the  scale 
parameter  3  and  shape  parameter  ex.  Using  these  results  we  demonstrated 
that  the  procedure  of  analyzing  the  effect  of  a  non-uniform  state  of  stress 
obtained  from  the  finite  element  method  can  be  post  processed  by 
statistical  reliability  formulation  accounted  for  size  effect,  resulting  on 
the  prediction  of  the  reliability  of  the  structure. 
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Figure  4.7  Element  Numbers  of  Half  Specimen  Model  at  Right 
Boundary  Region 
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Figure  4.8  Stress  Distribution  Values  of  Half  Specimen  Model  Based 
on  4-Point  Integration,  where  u'  =  0.0i",  at  Right 
Boundary  Region 


ELEMENT  *0.  1 

•U 

?  [  r,  7 

SIt.1 

SIC.» 

S1T.7  2 

41  = 

I . JU-5O0  1 

'|.ii24*750 

1 . 6*'4?t  2 

-i1 .  3  36?f -114 

0.21*57-03  -0.2320E-03 

4  JO 

1  .  "9  7S'l  0  i 

0.0  2  1*750 

1 . 0  O  Ml  C  —  *4  2 

-0 . 43646-04 

0.2uo7?-03  -0 . 3336E-03 

4?1 

4.  0  O65-!0.. 

1 . 0  *  1*750 

0.605  IF -0  2 

-0 . 7*  3oE -i'5 

0 . 2325E -0  3  — 0. 1M?=E— 04 

4?? 

6 .  .io=5:>o» 

0.  .120*750 

2. 6*o3t “02 

-0.»037E-'JS 

U.2291E-03  -0.57 4 6 £—04 

U?3 

0.0045000 

■1.02-1*750 

0 .60  *0E -42 

-0.111 u  ?  — 04 

0.2220E-03  -0 . 1 u  39£-0  3 

•4J4 

4. 09=5009 

0.»2<4»75o 

0.o335F-ii2 

-0 . 1436E— 04 

0.21 07E-03  -0 . 1 623E-03 

425 

il..  1035000 

0. o2'i*750 

'J.5055E-02 

-0 . 1 HS7E-04 

0.19U4E-03  -0 .237*E-o3 

4?6 

0  .  .100500  • 

0 .  (i2'i“  75H 

4 . 662*?-’)? 

-0.1  6*  7E-05 

0 .P308E-0 3  -".O534E-05' 

427 

n .  10  05  00 11 

0.02 4*750 

0.67205-02 

-0. 1  7.67E-05 

0.2271E-O3  -0 .20= 1 £-04 

u?f> 

0  .  r'005ftOu 

0.02-16750 

0.64*»f-42 

-0.251 5E-05 

0.21 *=? -03  -0. 3737E-04 

42<3 

n _  0O450OO 

■.'.<1246750 

0 . o  t  02E -0  2 

-0.3O70E-O5 

0.205bE-n3  -O.5H03E-O4 

430 

J  .  IWO5OO0 

11. C2  4 «75o 

0 . 55  I  5? -02 

-0. 3052E-05 

O.IS54E-03  -0.6726c -04 

431 

I. 124*751 

0 . 02**  75u 

o.»-i  o«e  -02 

-0 .  *3?3E— 14 

0 . ?u22£-o 3  -0 . 1 0 32E-02 

u>2 

n .  ■124  0  750 

0.02*6750 

1 . -'2a2E -  '2 

-0.51 65c -04 

0 .261 4f-03  -0.11 247.02 

433 

0.  .124*750 

0.024*750 

0.*5l *E-42 

-".  3277E-05 

0.2P*5f-03  -0. 12372-02 

434 

.  0  2  4  ?  7  5  y 

O.il2'i«75o 

4.07625-02 

0.56?5E-l)4 

0.321  IF -03  -U.13P3E-02 

4  35 

0.  124*750 

0 . 024*750 

0.0(,?oF-02 

0.2O59E-03 

0 .3  =  04? -0  3  -0.  1630E-02 

436 

0.  >29*75" 

0. 32u*75'J 

4.4UOI £-02 

0.2251E-03 

0.41  45F  -03  -.4.2108E-02 

437 

2.!' 24*75  0 

0.. 124*750 

).  1  1  *2E-  11 

0.  *»  3=E-03 

0 . 7b 1 4£  —0  3  -G.28"5c-02 

43* 

J  •  ‘2  4  *7  50 

0. 024*750 

4.  !l6?6*-02 

-0.  *55u£-05 

0.16*2?-o3  -0.97007-03 

430 

ij. 0243750 

it  .  ')  2  4  6  75  0 

3.457=7-42 

0.3261 c-o« 

G.16=2F-03  -1 .10307-02 

440 

0.024*750 

0.024*750 

j .  4  1 1. 3E - 42 

0 . 53o4c-04 

0.1 776? -03  -0.11007-02 

44| 

1.  >24*750 

0 . 0  2  4  *  7  5  0 

1. 3u2*c-o? 

C  .  1 761 : -03 

O.I*0Oc-O3  -0.11 7u?-02 

4  42 

0.024*750 

3. "24*750 

1.2377F— >2 

0. 2666E-0S 

0. 1 *“6£-03  -0.11 3 1 E -02 

443 

3.0246751) 

0.024*750 

C.2«*?E-43 

0 . 5 1 29E-03 

0.21  4=E  -1)  3  -0  .  “56 1  E-n 3 

4  44 

0. 024*75.1 

0.024*750 

0.2567E-)2 

0.16?6C-ll4 

0 . 94  7 1  £ -04  -0 .61 79?  —03 

'14  5 

0.024*750 

0.024*750 

4.21 1 2E-02 

O.5034E-04 

0. 9217E-04  — 0 ,0942c— 03 

4  46 

(1.  324*750 

0. "24*750 

1 . 1 5*o6“n2 

0.511 72-04 

0.9034E-04  -0.53B5E- 03 

447 

3.024*750 

0.0 24 *7 5" 

,  o-.^O? -4  3 

0. 1 439c-o J 

0 . 9  I  2 1 ?-04  -0.43*97.0 3 

44? 

0. 124*750 

0 . 0  ?  4*  750 

0.4 60**— 03 

0.  1  4  37E-03 

0.7359E-04  -0.3152E-03 

4  4  0 

0.024*750 

C..i24*  750 

4.=»1  75E  —  3 

o.*3?5E-')5 

■1. 3052£-04  -".20207-03 

430 

3.02**750 

o.r ?4"?50 

0.5372? -0  3 

".13  3  4  £  —  0  4 

".2364F-04  -0.17932-03 

'151 

0. 324*75'! 

'). 024*750 

1 .23"5?-"3 

0.6060?— 04 

4.  1.  =  *9E— 04  -4  .  1  358E-03 

432 

3.  '2  '1*7511 

l.i,2. 4  4  750 

1  .  S*  7 6?  - 0 4 

-.  3077c- 0  4 

0. 1  7277-4  4  -0.6723E-04 

453 

•0.  124*73'! 

•1.0  24*750 

0. 1  1  *  2E -0 1 

':.‘"i3«?-(,  3 

0.761  1  * -" 3 

||. ?*o  5c- 0  2 

4  ?4 

'■  .  724“  75,1 

0.1.  2  4  *7  50 

■>  .  44>.  J?  -02 

0.2? 40?-. 03 

4.4144E-03 

0.21 "77-02 

455 

0 .  12'!?  750 

0 .024*750 

..O.'lOC.-p 

"  .  ?r,5  7:-n  j 

". 3903c -03 

0.1630E-02 

4  5o 

1.. '24?75ii 

0."  2H-75II 

i.  .  ?  76  1 1  -"2 

0 . 60?  5E-04 

0.321  IE -0  3 

o.l 3* 2c -02 

'157 

3.024675) 

•j.o?U?7SG 

4. “51 35-  '? 

-0.326=7-05 

•l  .2005E-O  3 

0. 1237E-02 

45? 

...  .120*730 

I.. 1.24*750 

1  .•  2'.2£-o2 

.  5  1  5*7.04 

4.26I47-.J3 

i..l  1  2  3  c  -  0  2 

410 

■.  124*75  . 

•  . oi j-75  . 

.  .  !  '|C  -  6 

-0 . *i 1 4?-04 

0. 24232-03 

o.l 0  32E-0  2 

rnoi 

!0. 02*87501 

!0. 02<4875O| 

0. 24)736 -OSl 

'  .  5  1  2  a*  -i)  3 

".21477-U3 

f .955oc-03 

4  6  1 

.  2  i-./3'i 

o .  1!  2  4  4  7  3  .j 

...23/  7r  -  "2 

o.?6o4£  —  0  3 

0 . 10*6? -0  3 

0.1 1317-02 

4b2 

0 .024*750 

4.024*75  >4 

' . o  i2*E -0? 

0. 1  7S| 7- 03 

0. 1 H7  9? -0  3 

0.11 7  4  £ -0  2 

4  6  3 

:>.!>?  4*  75'> 

0.024*750 

0  .  41  *  2?  —  •"  2 

f‘.5364c  —  "4 

0.1  77b?  -(1  3 

0.1) 00c -02 

464 

32  3“  75(1 

0  .  .42  4  *  7  50 

1. J67=7-V 

0. 32657-04 

0  .  1  6  9  2  E  -  0  3 

0. 1 O30E-O2 

465 

’.024* 750 

0  .  >1 2  4  ?  7  5  0 

".  I3252-4 J 

-0. *7  =  47—15 

0 . 1 b«?£ -03 

fi.  =  7047-03 

466 

0. 02 ')  *75.1 

II  .  ')2  4*  7  50 

' .  '>6n2r  -  )  3 

0.14357-43 

0.73097—116 

0.31 40J-O3 

46? 

0.024* 750 

'I  .  <j  ?  4  -  7  5 11 

11.  °*  12* -"3 

0. ! -  3 * -  —  03 

0.911 b£-OJ 

0.43*5E-03 

4  6? 

0. ’24*750 

U.O?4*750 

(  .  15-5?— 42 

0  .  a  1  177-04 

0 .9060?  —(14 

u.53“0F-03 

460 

3. 12 o 375* 

• .020*750 

0.21 1 25-" 2 

O. 50 37E-04 

0.921 4E-0U 

0.5937E-03 

4  7  0 

1  .  .0  2  J  *  7  5  n 

II.OJU!'  75.1 

1 .2*60* -42 

4.1*2=7.04 

(1 . 04697  -04 

</.ol  75?-o3 

4  7  1 

..  •?  75,i 

0  .  1  *  4  <•  7  5  1 

0.S1-37E— •« 

f  .  707=7-114 

0.1 7246-04 

■|  .  b 7  0  3 £  -  0 4 

472 

■|.'l?7e  75" 

(1 .0244750 

>.??97?-P3 

n . 70007-04 

"  .  1  9*7t-''>4 

0.1 35bc-0  3 

473 

".  0  2  4  *  7  5  0 

fi. 02 '1-750 

4.5  36  0*  --'3 

0 . 1 3347-04 

9.2360C-O4 

0 . 1 7915.03 

4  74 

0. ,24*75  1 

'..<>2'i-7S3 

0.0)50^-03 

".*31 57-05 

0  .  5  4  4  7  E -04 

O.Jill  0E-O3 

475 

!•  •  .  1  2  4  ?  7  5  1 

".4*4-75  ' 

1.307) 5-4) 

r  .?ii*  =  ;-n? 

". 15237-02 

0.31  4*7.112 

176 

'•.  .2  1*7  5  • 

•1  .  3*4  7  5" 

-  '.-••ME-  ? 

•'.  1-  127-M 

-•.1.2734,  -  0  4 

1 . 1 2*47—03 

477 

■ .  2  4-/3  i 

'•■434-  76" 

■'.••2  3 

o.l  11  *:-"3 

•1  .  *'ISHf  -"  9 

p. 73777-05 

474 

7J-7  v,. 

0  .  1 3  1  -  75  ■ 

-0.1  .|*>2‘  -'•  3 

“  .  66  4  1  r  — I  4 

i.5b2iif-"5  -  ’.?*35c-‘l4 

Figure  4.9  Stress  Distribution  Values  of  Half  Specimen  Model  Based 
ON  1 -Point  Integration,  where  u'  =  0.001",  at  Notch 


ELEMENT  Hfl.  L  W  MilY  5  X &Z  SIC.*  5IT.Y2 
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Figure  4.10  Stress  Distribution  Values  of  Half  Specimen  Model  Based 
on  I -Point  Integration,  where  u' =  0.004",  at  Notch 


ELEMENT  NO.  L  *  51  f.T  $1&Z  SIC.»  Sir.YZ 

179  •.  t'>9S.)  ti1  :>.■.!•> /•nil'll  -'I.  I  J22-.-H4  •>.<'/ lit-UA  n.3ln/t-U9 

I  ^0  loas'-tift  O.oSTPE-o?  -0 . 1 3  I 88-03  0.27208-03  0.52835-Ob 

1*1  r,.'><»35i.n;)  ft.0S7'iftui>  7.°l  ObF -'>2  -O.llbbe— >3  0.2b2lt-03  0.an9bE-0a 

1*2  i.0995(oo  .->.ns7.)OP--i  ).'«J7<»E-')2  -0. 12*75-03  0.2o*>5F-0i  0.222bE-0« 

1 «  3  (i .  0  S  7  rif>  fi  ii  j.-'^'iSE-rt?  -0  .  t  3375-0  3  3.2b*«E-03  0.87t>b£-05 

11«  193510, ■»  ft .  057  mftO  0.95585-92  -n.13545-03  0.26995-03  0.29495-05 

MS  ". .'935000  (..OS7<;0"9  C.957o£-02  -0. 13585-03  0.27135-03  0.bba3£-0b 

1*6  ’.•.i'»95„un  r,. 1571001  ".on7if-r.g  “0.1)3S£“Q3  ft.2o21E-03  o.lbISE-oa 

T*7  r .  . . .  ■  • .  •*.« TTTTTTTo  ~ " >58* -Q.12738-Q3  1. 2*636-03  0.7896E-05 

1  -  *  .  .  .  ■  j  ",  5 . , , ,  "  "  i.  8  7  >.  ii  ft  o  ~ .  j**  “ft  .  1  337E-03  Q.?f»89f— 03  0.32915—05 

T5u~  "I .  .-aas./fti1  ; .  • .  5  7  • ,  ft  o  o  ■•i.9«.S.i»“ft>  -0 . 1 3825-03  0.2b04E-ns  0. 1  1565-05 

1  au  ",  -oaViiMi  i.-.ftS?  MiiiQ  7S>“|>2  -0 . 1  370E“n  3  0. 20985-03  0.27bbE-0b 

191  -  ti.’jVVVMXi  (,.0298750  j.<52Sut“ii2  -C.133IE-03  0.2oft«E-03  0.1649E-03 

192  O.J99500'  <>.  9?a*759  <,.90378-02  -0.1P15E-03  O.2578F-03  0.19755-03 

193  0.'J®35CiOO  o.129«75‘»  ft .  8*9 1  f -12  “ft.H05E“03  0.2572E-03  0.1  133E-O3 

I  au  ft.u99SC0u  0.029*750  f:.H«01F-O2  -O. 10255-OS  O.2579F-03  (i.7()aQE“0« 

19S  .1.0995010  ft  .  0  2  9  *  7  5  0  ft  .  ^7^ap.ftg  -0.98305-00  0.25765-03  o.23(J6E-0u 

nr  v  .  ‘ 1  u  3  S  - 1  .J  o  <i .  ii  2  a  a  7  S  0  ft.9G008.i2  -9.12*55-03  0.2S37E-03  0.231  3F-03 

197  j. 7995101  0.0208750  ft .  m  7  66?  “02  -O.I121S-&3  ». 25295-03  (I.192SE-03 

198  ft.  ftftaso-fti  ft. ft?: 1875')  ft. 75-92  — ft  .  9  8  7  6  E — ft  a  ft  . 2S?bF  -0  3  0.19185-03 

1  a<9  ft.ft9a5.10-1  0 . 0 2 9 8 7 5 0  9.65305-12  -<> . *9705-04  U.2533E-03  0.8bl 75-00 

?oo  ft.  lOQS'ftft  9,ft?4*750  -1 .  a  ft  a  |  C  -  ft  ?  -o.«5l95-oa  0.2538E-O3  0.2hbbE-0a 

201  ;...  >995ftjj  y.  029-751'  ft  ,-.M5!“i'2  -0.11595-03  1.2UU9E-03  0.29995-03 

202  0. r  99500- 1  u.w?9»7S  )  ■..-•3v3r—ft?  -ii.975bE-oa  i>.2aS5E-«3  0.23U8E-03 
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Figure  4.11  Stress  Distribution  Values  of  Half  Specimen  Model  Based 
on  1 -Point  Integration,  where  u'  =  0.001",  at  Right 
Boundary  Region 
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Figure  4.12  Stress  Distribution  Values  of  Half  Specimen  Model  Based 
on  1 -Point  Integration,  where  u'  =  0.004",  at  Right 
Boundary  Region 
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0.009071 
0.009358 
0.009496 
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0.009750 
0.008758 
0.008491 
0.008179 
0.007840 
0.007500 
0.007195 
0.006968 
0.006855 


0.008717 


U'  =  0.004' 

U'  =  o.or 

0.036280 
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0.037420 

0.090350 

0.037970 

0.094830 

0.038190 

0.095370 

0.038280 
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0.035030 
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0.033960 

0.084860 
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0.031350 

0.078260 

0.029980 

0.074810 

0.028760 

0.071720 

0.027840 

0.069400 

0.027880 

0.068220 

0.034850 

0.087045 

TABLE  4.  2 


AVERAGE  STRESS  FOR  FULL  SPECIMEN  MODEL  WITH  NOTCH 


ELEMENT 

tfy  ,  Msi. 

NUMBER 

WIDTH,  IN. 

U'  =  o.or 

163.  166 

0.095000 

0.091560 

164,  167 

0.095000 

0.094940 

165.  168 

0.095000 

0.095670 

171,  465 

0.024875 

0.087600 

174,  468 

0.024875 

0.084860 

177,  471 

0.024875 

0.081660 

180,  474 

0.024875 

0.078170 

183,  477 

0.024875 

0.074670 

186,  480 

0.024875 

0.071540 

189,  483 

0.024875 

0.069200 

192,  486 

0.024875 

0.068030 

PBAR* 

0.078025 

0.078025 


■tw  |l. 


TABLE  4. 3- A 

STATISTICAL  DATA  SUMMARY 
FOR  PLATE  WITH  NOTCH ,  «  =  I0,  jl  =  0.1,  kc  =  0.001 


MODEL 

HALF 

FULL 

u#=o.oor 

U -0.004" 

U'=0.0I" 

U’=0.0I" 

PBAR* 

0.008717 

0.034850 

0.087045 

0.087025 

h 

0.030582 

0.030749 

0.03II10 

0.028880 

FOR 

PLATE  WITH  N 

OTCH ,  <x  =  2C 

K  fi  =  0.1,  kc 

=  0.0001 

MODEL 

HALF 

FULL 

u*=o.oor 

U'=0.004" 

U'=o.or 

U'=0.01" 

pbar* 

0.008717 

0.034850 

0.087045 

0.087025 

h 

0.029122 

0.029298 

0.029721 

0.028826 

TABLE  4.3-C 

STATISTICAL  DATA  SUMMARY 


FOR  PLATE  WITHOUT  NOTCH ,  «  =  10,  =  0.1,  kc  =  0.001 


TABLE  4.  3-D 

STATISTICAL  DATA  SUMMARY 
FOR  PLATE  WITHOUT  NOTCH ,  a  =  20,  fi  =  0.1,  kQ  =  0.000 f 


Average  boundary  stress 
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APPENDIX  A 


STATISTICAL  STRUCTURAL  RELIABILITY  PROGRAM 


c*»» . . . . . . . 

c: 

. . . . . . . 

c 

c.  . . . . . . 

PROGRAM:  RELIABILITY  CALCULATION  OF  PLATE  MODEL  BASED  ON 

LINEARIZE  SIZE  EFFECT,  NONUN l FORM  UNIAXIAL  STRESS. 
RELIABILITY  FOAMULA  T  ION  IS  EXPRESSED  IN  THE 
*  STANDARD  VIE  I  BULL  FORK  *. 

DATA  INPUT;  OUTPUT  OF  FINITE  ELEMENT  COOE  NAME  *  NIKE2D" 


VARIABLE:  KC  a  MATERIAL  CONSTANT  FOR  PLATE  THAT  IS 

SUBJECTED  TO  UNIFORM  TENSILE  STRESS.  IN 
THIS  CASE.  ASSUME  IN  EACH  ELEMENT  THE 
LOCAL  TENSILE  STRESS  IS  UNIFORM. 

BETA- a  SCALE  PARAMETER  FROM  EXPERIMENTAL  RESULT 
FROM  SEVERAL  NUMBERS  Of  SAMPLE. 

ALPHA  a  SHAPE  PARAMETER  FROM  EXPERIMENTAL  RESULT 
FROM  SEVERAL  NUMBERS  OF  Sample. 

PBAR  a  AVERAGE  STRESS  ALONG  BOUDAPY  side  of  the 
PLATE  THAT  SUBJECTED  TO  THE  EXTERNAL  FORCES. 

IT  MUST  BE  CALCULATED  SEPARATELY  AFTER  obtaining 
the  STRESS  RESULTS  FROM  "NIKE2D*  PROGRAM. 

NE  a  NUMBER  OF  ELEMENTS  OBTAINING  FROM  MESH  GENERATOR 
PROGRAM  NAME  "MAZE”  WHICH  IS  THE  PROGRAM  THAT 
CREATE  AN  INPUT  FILE  TO  BE  USED  IN  "NIKE20" 
PROGRAM.  ALSO  BOUNDARY  ELEMENT  NUMBERS  CAN  BE 
CREATED  FROM  "MAZE"  PROGRAM. 


DOUBLE  PRECISION  F ,F4 ,FF,P( »000>, RA( 1 OOO) ,BETAA 
DIMENSION  XL^23G8'  ,XWC236B)  ,SI&tL23(jB> 

ReAL  KC 

byte  aisd; 


onoonoooo 


type*. 'enter  pile  containing  stresses  data  inrut* 

READ(5, 5}  A 
5  F0RMAU3OA1) 

OPEN(UNITsA.  FILE*A.  STATUS* ’OLD *) 

OPENTUNIT**,  FIL£*'SS.DAT*,  STATUS* 'NEW ) 

WRITE  I A,  1 ) 

1  FORMAT (3x, 'DATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT 
♦UNIAXIAL  TENSILE  STRESS') 

TYPE*, 'ENTER  1  FOR  PLATE  WITH  NOTCH 
»  OTHERWISE  FOP  PLATE  WITHOUT  NOTCH* 

READ*.NC 

IFfRC.EQ.l)  THEN 

WRITE (6 , 2) 

ELSE 

«RITE(6,22> 

END  IF 

8  FORMAT (20* , *  FOR  PLATE  WITH  NOTCH*) 

22  FORMAT (1  a*. 'FOR  PLATE  WITHOUT  NOTCH*) 

TYPE*, 'ENTER  THE  VALUE  OF  DISPLACEMENT  BOUNDARY  CONDITION' 

REAO* ,U 
NRlTE(fc.3)  U 

TYPE*. 'ENTER  1  FOR  ENURE  SPECIMEN  INPUT 
»  OTHERWISE  FOR  HALF  SPECIMEN  INPUT* 

READ*. NS 
IF (NS.ES.l )  THEN 
WR  IF  £  16.23) 

ELSE 

WRITE (b, 24) 

END  IF 

23  FORMAT (I5X. 'ENTIRE  SPECIMEN  MODEL  INPUT*) 

24  FORMAT* I6X. ’HALF  SPECIMEN  MODEL  INPUT') 

3  FORMAT (7X. 'WHERE  THE  DISPLACEMENT  BOUNQARY  CONDITION  IS  ',F8.«) 
TYPE*. 'ENTER  The  KC  VALUE' 

READ* ,KC 

TYPE* , ' ENTER  THE  SETA  VALUE' 

PEAO*,BETA 

TYPE*. 'ENTER  THE  ALPHA  VALUE' 

REAO*. ALPHA 
AlPHA2*AlPMA-»2, 0 
FFaO.O 

type*. -Enter  the  pbar  value* 

READ»,PBAR 

type*. 'ENTER  THE  NUM8FR  OF  ELEMENTS' 

PEAD».N£ 

REAO  ELEMENT  LENGTH,  WIDTH,  ANO  TENSTLE  STRESS  FROM  INPUT  FILE 

calculate  setaa  which  is  beta  that  changed  with  geometry  of  the 

SPECIMEN  MODEL  ANO  CERTAIN  LOADING  BOUNDARY  CONDITION  BUT  IT  IS 
INDEPENDENT  OF  THE  VALUE  OF  LOAD  UNDER  Th£  SAME  BOUNDARY  CONDITION 
00  200  1*1, NE 

REAO  14,  20  )  XU.r>.X»U>,SI<irCI) 

JFtSrGYID.GT.O.O)  THEN 
F-(XL( I)*XWf I) /KC  *  *2.0) 

»  *(31GT<I)«*»ALPHA2> 

FFsfPtF 


ci  n  o 


END  IF 

300  CONTINUE 

IF  f  NS . EQ . 1 )  THEN 
Ff  »FF 
ELSE 

Ff sff*2.0 
END  If 

FA*FF **<- I ./ALPHAS )»P8AR 

B£TAA’BETA»»(ALPHA/ALPMA2)*PA 

NRITEId.lO)  KC 
tO  F0RMATT8X, 'KC  a',Ft3.S) 

NRITEC6.lt)  SETA 
tl  FORMA T(6X. 'SETA  a'.F13.S) 

NR  l TE(6« 12)  ALPHA 
IS  FOR-AKSX, 'ALPHA  *'»f  ti.8) 

NRirE(6,13)  PBAR 
15  FORMAT (6X.'PBAR  *',F13.at 
NR  I  T  E  (  6 ,  1  4 1  HE 
t»  FOR«»r(»X. 'HE  a*. 14) 

WRITER, 15)  ALPHAS 
t5  FOR^A  f  (4*  ,  'ALPHAS  J'.fU.J) 

NR  I TE( 6,25 )  SETA A 

TYPE*. 'ENTER  THE  FIRST  EXTERNAL  LOAD* 

REA0»,P(1) 

PtNCsPf 1 )/30. 

RELrASILITy  CALCULATION  BASED  ON  LINEARTZE  SIZE  EFFECT. 
DO  100  J*1,100 

RACJ)=OEXPL-<Pf J1/#ETAAI« ‘ALPHAS) 

wpire(6.3oi  j.p(j).RAU) 

PUnt»PO>vPlMC 
100  CONTINUE 

SO  FORMAT («X.2F11.0,Et2.0) 

3S  FoRhAT'SX. 'BfTAA  =  ' ,Fi3.9) 

30  FORMAT'S*. 'X  s', 14. 3*. 'P  *' .F13.8.3X, 'RA  *'.Fl3.a) 

STOP 


APPENDIX  B 


DATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 

data  output  non  uiNfiAurze  Size  effect  uniaxial  tensile  stress 
for  Plate  with  notch 

where  t«e  displacement  boonoart  condition  is  o.ooio 


half  specimen  model  input 


KC 

s 

0. 

90010000 

SETA 

s 

0. 

10040000 

Alpha 

s 

20. 

00000000 

P8AR 

= 

0. 

00871700 

NC 

s 

482 

ALFHA2 

s 

22. 

ooooaooo 

SETAA 

= 

0. 

02912212 

J 

s 

t 

P  3 

0.00871700 

RA  3 

1 .00000O00 

J 

3 

2 

P  3 

0.00900757 

RA  3 

1.00000000 

.1 

= 

I 

P  3 

0.00924513 

RA  s 

1  .00000000 

J 

3 

'1 

P  3 

0.004584  7 y 

Ra  3 

t  .0000001)0 

J 

3 

5 

P  3 

0.00987927 

RA  3 

1.00000000 

J 

3 

(j 

P  3 

0. O10tb9«3 

RA  3 

1.00000000 

J 

: 

7 

P  3 

0.01046040 

RA  3 

1 .00000000 

J 

= 

A 

P  3 

0.01075097 

RA  3 

1.00000000 

J 

3 

9 

P  3 

0.01 1 >4153 

RA  3 

1.00000000 

J 

3 

to 

P  3 

0.0  1  1  332  H* 

RA  3 

1 .90000000 

J 

5 

1  t 

P  3 

0.01 162267 

RA  3 

1.00000000 

J 

3 

12 

P  3 

0.01191323 

RA  3 

1 .03000000 

J 

3 

1  3 

P  3 

0.01220580 

RA  3 

1.00000000 

J 

3 

14 

P  3 

0.01249437 

RA  3 

0.99999999 

J 

= 

15 

P  3 

0.01278493 

RA  3 

0.99999999 

J 

3 

1  b 

P  3 

0.01307550 

RA  3 

0.99999995 

.1 

3 

1  7 

P  3 

0.01536607 

RA  3 

0.99999996 

J 

3 

1H 

P  3 

0.01 365663 

RA  3 

0.99999994 

t 

3 

19 

»  3 

0.01 394720 

RA  3 

0.99999991 

I 

3 

20 

P  3 

0.01423777 

RA  3 

0 . 99999485 

J 

= 

21 

P  3 

0.01452533 

RA  3 

0.99999977 

J 

3 

22 

P  3 

0.01481890 

RA  3 

0.9999906$ 

J 

3 

2  3 

H  3 

0.01510947 

RA  3 

0 • 99999446 

J 

3 

2'J 

P  3 

0.01540 Oo j 

RA  3 

0.99949918 

J 

5 

as 

P  3 

0.01569060 

RA  s 

0.9  9499877 

J 

3 

2D 

y  - 

0.015981 1 7 

RA  3 

0.49999815 

J 

= 

2T 

o  s 

0.0162/1 73 

RA  s 

0.99909725 

] 

3 

28 

P  3 

0.0 1 6  5  o  2  3  0 

RA  s 

0.99999595 

J 

3 

29 

P  3 

0.916852*7 

RA  s 

0 . 94999406 

J 

3 

30 

0  3 

9.0! 714343 

9  A  3 

6.99999135 

J 

3 

31 

P  3 

0.0  17  '1  3 400 

9  4  3 

0.99995747 

J 

3 

13 

P  3 

0.01  772  157 

RA  s 

0.99998198 

J 

3 

3 1 

P  3 

0.91801JI3 

RA  s 

0.99907423 

J 

3 

34 

3  - 

0.01830570 

RA  s 

0.99996336 

J 

3 

45 

P  3 

0.01554627 

RA  3 

0.99994819 

J 

3 

3b 

P  3 

0.0! 858643 

RA  3 

0.99942712 

J 

3 

37 

P  3 

O.0191 7790 

RA  s 

0.99489804 

J 

3 

38 

P  3 

0 . 9 1  9  4  a  7  9  7 

RA  s 

0.999*5806 

J 

3 

39 

P  3 

0.01975853 

RA  s 

0.99980337 

J 

= 

49 

P  3 

0.02004910 

4  A  s 

0.99972890 

J 

3 

41 

P  = 

0.  >21133467 

RA  s 

0.99962797 

J 

3 

9.3 

°  3 

0.9296  5  125 

RA  s 

0.99949175 

J 

3 

«5 

P  3 

0.9 9 9  42  9 8(1 

RA  3 

0.99930870 

J 

3 

44 

P  3 

0.02121 1 57 

»A  3 

0.99906373 

J 

= 

4S 

P  3 

0.02150195 

RA  : 

0.99673722 

J 

3 

46 

J  - 

0.991  7925  > 

RA  r 

0.99830378 

J 

3 

47 

P  3 

0. O229m$07 

RA  3 

0.99773061 

.1 

3 

4« 

P  3 

o.«22375o? 

RA  s 

0.99007S58 

J 

3 

44 

O  - 

1 .02266420 

RA  s 

0 . 49598(1  76 

64 


0  T*  OUTPUT  fop  Lt*£A *IZE  Size  EFFECT  UNIA*IAC  tensile  stress 
FOR  PLATE  *rTH  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  XS  O.OIOO 

„  HALF  specimen  mooel  input 

KC  3  0.00010000 

BETA  3  0.10000000 


Alpha 

« 

20 

.oooooooo 

PBA* 

3 

0 

.08704900 

HE 

3 

*82 

ALPHAS 

3 

22 

.00000000 

ae  taa 

s 

0 

.02972093 

j 

= 

1 

P  3 

0..  10671  700 

9  A  3 

1 .OOOOOOOO 

j 

S 

2 

P  3 

0.:)  .1900757 

9  4  3 

1.00000000 

j 

2 

J 

P  3 

0.094304  |  3 

94  s 

1  .oooooooo 

j 

s 

a 

P  3 

9.094SAS70 

94  3 

1.00000000 

j 

2 

5 

P  3 

1,00907927 

94  3 

l  .OOOOOOOO 

j 

3 

o 

P  3 

0.0101o98 J 

94  3 

1.00000000 

j 

= 

7 

P  3 

0. 01096040 

94  s 

1 .OOOOOOOO 

j 

- 

9 

P  3 

0.01075 9 °7 

9  A  s 

1 .OOOOOOOO 

j 

3 

9 

P  3 

0.01  1  OU|<53 

44  3 

t. OOOOOOOO 

.i 

5 

to 

P  3 

1.01 153210 

9  A  3 

1 .oooooooo 

j 

S 

1 1 

P  3 

0.01 1 42267 

94  s 

1.00000000 

j 

- 

12 

P  3 

O.o 1 1 9 1  523 

4  A  3 

1.00000000 

j 

= 

1  s 

P  3 

0.0 1 220  390 

9A  3 

1  .oooooooo 

j 

S 

1  <1 

P  3 

0.01209937 

9  A  s 

0.99999999 

j 

15 

D  S 

0.01279493 

94  a 

0.99999999 

j 

s 

to 

P  3 

0.01 307550 

94  3 

0.94999999 

j 

3 

17 

P  3 

0.01 336607 

94  3 

0.99999998 

j 

3 

19 

P  3 

.0.01  365663 

94  3 

0.99999996 

1 

3 

19 

P  3* 

0.U1  3997<>i> 

94  3 

0.9999999a 

j 

3 

JO 

P  3 

0.0t9<>3777 

94  s 

0.99999991 

j 

3 

Jt 

6  z 

0. 01052433 

44  3 

0.99999905 

j 

3 

22 

P  3 

0.01441 490 

94  3 

0.99999978 

j 

3 

23 

P  3 

0.01510907 

44  3 

0.999999o6 

j 

S 

J« 

P  3 

0. 01500003 

94  3 

0, 9399994a 

i 

3 

25 

P  3 

0.01569060 

94  3 

0.99999921 

j 

S 

26 

3  z 

0.015981 1 7 

44  3 

0.99999002 

j 

2 

27 

P  3 

0.016271 73 

9  4  3 

0.99909895 

j 

20 

P  3 

O.o  1  o5o J  So 

94  3 

0.9999Q7UI 

j 

J9 

P  3 

O.u16°5?07 

4  A  s 

0.99999620 

j 

2 

JO 

P  3 

O.oi 710343 

94  3 

0.99999a 9 7 

j 

2 

31 

P  3 

0.01793900 

94  s 

0.99999200 

j 

S 

3? 

P  3 

0.01 772457 

9  A  3 

0.99990090 

j 

2 

S3 

P  3 

0.01001513 

94  3 

0.99990353 

j 

2 

P  3 

0.01930570 

94  s 

0.99997658 

j 

2 

35 

P  3 

0.01459627 

44  s 

0.99996600 

j 

2 

3b 

P  3 

0.01444643 

94  s 

0.90905393 

j 

3 

37 

P  3 

0.01917740 

94  3 

0.9999390a 

j 

3 

38 

P  3 

0.0  19-10  797 

94  3 

0.93990029 

j 

3 

39 

P  3 

0.OI97535J 

94  s 

0.99987939 

j 

2 

«o 

P  3 

0 . 0  2  0  9  U  9 1  o 

94  3 

0.99902675 

j 

2 

HI 

P  3 

0. 02053967 

94  3 

0.99976229 

j 

2 

02 

P  3 

0.02063023 

94  3 

0.90967518 

j 

2 

?  3 

1. 09032000 

44  3 

0.99955017 

j 

2 

u« 

P  3 

0.02121 1 37 

94  3 

0.93990150 

j 

2 

45 

P  3 

0.02150193 

94  s 

0.9991 92Ho 

j 

2 

4b 

P  3 

a.o2 i 79450 

94  3 

0. 99091570 

j 

2 

47 

P  3 

0. 0920430  7 

94  3 

0 . 930599 1 b 

j 

2 

46 

P  3 

0.09?373o5 

94  3 

0.99006620 

•i 

2 

49 

P  3 

0 . Q  22  60  42  0 

U  4  3 

0.997-13222 

data  OUTPUT  FOR  LINEARIZE  SIZE  effect  uniaxial  tensile  stress 
FOR  PLATE  WITH  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  O.OONO 

half  specimen  model  input 


xc 

S 

0. 

00050000 

seta 

3 

0. 

10000000 

ALPHA 

s 

20. 

00000000 

P6AR 

s 

0. 

03A85000 

NE 

3 

462 

ALPHAS 

3 

22. 

00000000 

8ETAA 

3 

0. 

02929771 

J 

3 

1 

P  s 

J.  00871  7(10 

PA  s 

1.00000000 

J 

3 

2 

P  s 

0.40900737 

PA  3 

1 .0000(1000 

J 

3 

J 

P  3 

0. J092931  5 

PA  3 

1 .00000000 

J 

3 

a 

P  2 

0.0095  5570 

PA  3 

1 .00000000 

J 

3 

S 

P  2 

0.00947927 

PA  3 

1.00000000 

J 

3 

6 

P  2 

0.91018985 

PA  3 

1 .0001)0000 

J 

3 

7 

p  • 

0.11046040 

PA  3 

1 .0001)0000 

J 

3 

8 

P  s 

0.0  to  7509 7 

PA  3 

1 .000(10000 

J 

3 

* 

P  3 

0.01114155 

PA  3 

1.40000000 

J 

3 

1  n 

P  s 

0.01153210 

PA  3 

1.00000000 

J 

3 

11 

P  s 

0.01 1622b? 

PA  3 

1.03000000 

J 

3 

12 

P  3 

0.0! 191525 

PA  3 

l .00000000 

J 

3 

1  I 

P  3 

0 . 0 122035" 

PA  3 

1 .00000000 

J 

3 

t 

P  3 

0.01299457 

PA  3 

0.99999999 

J 

3 

IS 

P  3 

0.01274495 

PA  3 

0.99999099 

J 

3 

1  6 

P  • 

0. It  307550 

PA  s 

0.9999099)} 

J 

3 

1  7 

^  3 

0.01336607 

PA  3 

0.99990997 

J 

V 

18 

P  3 

O.oi 345495 

PA  3 

0.99090999 

J 

3 

19 

P  3 

0.01  59.1720 

PA  3 

0 .99999902 

J 

3 

?a 

P  s 

0.  11425777 

PA  3 

0 . 0990004  7 

J 

3 

>i 

P  s 

0.01  152833 

PA  3 

0 • 999999H0 

J 

3 

?2 

P  3 

.).•■' 1  451  49» 

PA  3 

0.09000060 

J 

3 

2  5 

P  3 

0.0151 094 J 

PA  3 

0.09090055 

J 

3 

2  4 

r»  - 

0.015  10005 

Pa  s 

0.999O9»2« 

i 

3 

55 

P  3 

0.OI54O06.I 

PA  3 

0.99900402 

J 

3 

26 

P  3 

0.IH59H1  1  / 

PA  3 

0.94000456 

J 

3 

27 

P  3 

l.0lo?7i 73 

PA  3 

0. 09909759 

J 

3 

?“ 

P  3 

0.  11654250 

PA  3 

0.90909645 

.1 

3 

jo 

P  3 

0.  Jlb452a7 

PA  3 

0.99999479 

J 

3 

5  j 

P  3 

•1.01  7  11  54  5 

PA  3 

0.99999242 

J 

3 

51 

P  3 

0 . 0  1  7  4  5  1 1.1 0 

»  A  3 

0.99998005 

J 

3 

52 

°  3 

0.01  772457 

PA  s 

0.9099842! 

J 

3 

33 

P  3 

0.01501515 

PA  3 

0. 00007742 

J 

3 

54 

P  3 

0.01450570 

PA  3 

0.99906700 

J 

3 

35 

P  3 

0.0  1  454427 

PA  3 

0.06095460 

J 

3 

5b 

P  3 

0.01654645 

PA  3 

0.99095615 

.J 

3 

37 

P  3 

0.0191  7740 

PA  3 

0.99901067 

J 

3 

38 

P  3 

0.0  1  <J0r,797 

PA  3 

0.49047564 

J 

3 

no 

P  3 

0.01975555 

PA  s 

0.09042772 

J 

3 

'4  0 

P  3 

0.0200441O 

PA  3 

0. 90076244 

J 

3 

4  1 

P  3 

0. 02M344  7 

PA  3 

0.99947405 

J 

3 

4  2 

P  3 

0. 42463125 

PA  3 

0.99055" 7  0 

J 

3 

4  5 

P  3 

0.  )21920<(i 

PA  3 

0.99039451 

J 

3 

4  4 

P  3 

0.02121157 

PA  3 

0.0901 7067 

J 

3 

4  5 

P  3 

0.021 501 95 

PA  3 

0.49949557 

; 

3 

4* 

P  3 

0.021  7-9  25  0 

p  A  3 

0.9445! 575 

j 

3 

«  1 

P  3 

0  .  "22  1 4  30  7 

PA  3 

0.99401 1  1  7 

j 

3 

au 

P  S 

0.  122  57  563 

PA  z 

0.99751976 

j 

- 

»o 

t»  3 

•1.4  »2  61 14.. 

Pd  3 

0.996441 50 

OAT*  OOTPur  PO*  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILF  STRESS 
FOR  PLATE  VfITHOHT  NOTCH 

WHERE  THE  DISPLACEMENT  gOUNOART  CONDtnON  IS  O.OUIO 
ENTIRE  SPECIMEN  MODEL  INPUT 
KC  s  0.000(4040 
SETA  :  0.10044000 

ALPHA  3  24.00040000 

PSAR  3  0.009SIZ0O 

N£  3  200 

ALPHAS  3  22.00400004 

SETAA  3  0.05180000 


J  3 

l 

P  3 

0.04*71 700 

P  A  3 

1 .00000000 

J  3 

2 

P  3 

0. •■10900757 

PA  3 

1.00000000 

J  3 

J 

P  3 

•I.  4692°*  1  J 

PA  3 

1 .00000000 

J  3 

4 

P  3 

0 .0095447" 

PA  3 

1 .00000000 

J  3 

S 

P  3 

0 • 0094 7927 

PA  3 

1.000 00000 

J  3 

0 

P  3 

O.OlUlo***? 

PA  3 

1  .000000(70 

J  3 

7 

P  3 

0.01 044040 

PA  3 

1.00000000 

J  3 

A 

f>  - 

0.01 0 75 09 7 

PA  3 

1.00000000 

J  3 

0 

P  3 

0.01 10415? 

PA  3 

l  .1)0000000 

J  3 

10 

P  3 

o.oi i ??2io 
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IS 
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[0 

P  3 
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PA  3 
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20 

P  3 
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J  3 

21 

P  3 
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J  3 

j; 
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PA  3 

1 .000000 00 
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3a 
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PA  r 

1  .1)0000000 

.1  3 

2S 

P  3 

0.01S6904( 

PA  3 
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1  3 

2b 
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J  3 

27 

O  - 

0.01o27l 75 

PA  3 

1  .0  9000000 

J  3 

28 

P  3 
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PA  3 
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J  3 

29 

P  3 
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1  3 

JO 

P  3 

0.  >1  7145-13 

PA  3 

1 .00040904 

J  3 

SI 

O  2 

O.'H  745400 

PA  3 

1 .  o  o  o  o  o  u  0  o 

I  3 

S2 

P  3 

O.ol 772457 

PA  3 

1 .00000900 

J  3 

ii 

P  3 

n.0160lj|5 

PA  3 

1 .09000000 

J  3 

J* 

P  3 

0.91050570 

PA  3 

1  . ooooooou 

J  3 

S5 

P  3 

0.01  45*427 

PA  3 

1 .00009000 

J  3 

So 

P  3 

). 116**635 

PA  3 

1.04000004 

J  3 

S7 

a  2 

0.0191 7740 

PA  3 

1  .O0O40OUO 

I  3 

S8 

P  3 

0. 01946747 

PA  3 

1  .0000090  0 

J  3 

SO 

U  2 

0.01 775653 

PA  3 

1 .04090000 

J  3 

40 

P  3 

o.  ;  P<>  0  4  9  ( •> 

PA  3 

1 .ooooooou 

J  3 

<U 

P  3 

0. (2055967 

PA  3 

1 .01000900 

J  3 

42 

D  - 

o .729,5123 

P  A  3 

1  .40100009 

J  3 

43 

O  » 

■1 .0  20  4?  7*9 

PA  3 

l  .04090400 

J  3 

44 

P  - 

0 . 9 2 1 2  l  1  3  7 

PA  3 

1  .  90000900 

J  3 

«5 

P  3 

0.021  3<»l  73 

PA  3 

1.00940000 

J  3 

46 

P  3 

0.021 792,0 

PA  3 

0.99999099 

J  3 

47 

P  3 

0.02*06  507 

PA  3 

0 . 99399°99 

1  S 

uS 

P  3 

0.0223756? 

PA  3 

|)  .  04000060 

J  3 

4.9 

U  2 

o.o 2266  12" 

PA  3 

o . 90O4OQPO 

OATt  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT 

oniaxml  tensile  stress 

FOR  PLATE  NITHOUT 

NOTCH 

WHERE 

THE 

DISPLACEMENT  BOONOABY  CONOITIOM  IS  0.0100 

ENTIRE  SPECIMEN  MODEL 

INPUT 

KC 

s 

0. 

OOOl 0000 

BETA 

s 

0. 

10000000 

ALPHA 

s 

20. 

ooooaooo 

PBAR 

s 

0. 

09091 000 

HE 

s 

200 

ALPHA  2 

s 

22. 

00000000 

8CTAA 

s 

0. 

05180009 

J 

s 

1 

P  3 

'1.00871  7yO 

PA 

3 

1.00000000 

J 

s 

2 

P  3 

0  .(1090  0  7S 7 

PA 

3 

1 .00000000 

J 

s 

5 

P  3 

0. -109294  |  J 

PA 

3 

1 .onnoooon 

J 

s 

‘l 

P  3 

0.9095887-7 

PA 

3 

1.00000000 

1 

s 

5 

p  3 

0.0098792/ 

RA 

3 

1.00000000 

J 

s 

b 

P  3 

0.01016999 

RA 

3 

1.00000000 

J 

s 

7 

P  3 

0.0109604(1 

Rd 

3 

1.00000000 

J 

s 

n 

P  3 

0.010/5397 

RA 

3 

1  .(>1000000 

J 

S 

9 

P  3 

0 . 0  l  1  0  4 1  5  5 

PA 

3 

1 .00000000 

J 

= 

P  3 

0.01155210 

RA 

3 

1  .oooooooo 

J 

s 

1  I 

P  3 

0.4 1 162267 

RA 

3 

l .no jnoooo 

J 

= 

12 

P  3 

0.01191525 

RA 

S 

1.00000000 

J 

S 

t  5 

P  3 

0 . 1 1 220  5 Ad 

RA 

3 

1 .oooooooo 

J 

s 

19 

P  3 

rt.11 299457 

RA 

3 

1 .oooooooo 

J 

3 

IS 

P  3 

0  .  ()  1  2  7  A  9  9  5 

RA 

3 

1.00000000 

J 

S 

tb 

P  3 

0.01  50  755-7 

RA 

3  . 

1.00000000 

J 

S 

t  7 

P  3 

0.01 5  5ob1 7 

RA 

3 

t.nooooourt 

J 

» 

1 1 

P  3 

3. 0  1  5.65665 

PA 

3 

1.00000000 

.1 

S 

19 

P  3 

O.ul  594  72,1 

RA 

3 

1 .07000000 

J 

s 

20 

P  3 

-1.01  92  5  777 

R  A 

3 

1 .05000000 

J 

S 

21 

P  3 

0.01952655 

RA 

3 

1 .oooooooo 

J 

s 

22 

P  3 

0.01 951490 

PA 

3 

1.00000000 

1 

3 

25 

P  3 

(1.1/(8(064  7 

PA 

3 

1 .00000(10(1 

J 

S 

2  '4 

P  3 

0. 01540115 

PA 

5 

1 .000(10000 

1 

3 

js 

P  3 

0.0  1569(1 ->n 

RA 

3 

1 .00010 non 

J 

3 

29 

P  3 

0.11S4AH  / 

PA 

3 

1  .0500001)0 

J 

3 

27 

P  3 

•1.0162/173 

RA 

3 

1 . 0  -1 0  0  0 1 0  0 

J 

3 

28 

P  3 

0  .  1  6  3b  2  5  i 

RA 

3 

1 .07000001 

.1 

3 

2° 

?  3 

').  01  6  852  57 

RA 

3 

1 .oooooooo 

J 

3 

in 

P  3 

7.01719545 

RA 

3 

1.3000000(1 

J 

s 

31 

P  3 

1  .  I  1  7-4  59  1  1 

P  A 

3 

1 . 0  0  0  0  0  0  0  0 

J 

3 

32 

M  S 

0.11  772  457 

PA 

5 

1  .ynonoo«o 

J 

3 

33 

P  3 

0. 01  Art  151 5 

RA 

3 

t  .01000000 

J 

3 

3« 

P  3 

0.0 1,8  5(1 87(i 

RA 

3 
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J 

3 

3S 
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RA 

3 
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= 

3b 
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0  A 

5 
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= 
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PA 

3 
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J 
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PA 

3 
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J 
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PA 
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J 

= 
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C  - 
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PA 

3 

1  .OOOOOOOO 

J 

3 

41 

P  3 

0.12055967 

RA 

3 

1.00000000 

J 

3 

92 
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PA 

3 

l .oooooooo 

J 

= 

43 

P  3 

0.  12092  18. 1 

RA 

3 

1 . 0)00 0000 

J 

3 

44 

P  3 

0 . 1 4  1  2  1  1  57 

RA 

s 
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J 

s 

<J5 

P  3 

I.  021  SOI  15 

PA 

3 

1.070-1000  0 

J 

3 

49 

0  • 

-1.191  79950 

RA 

s 

0 . 09909OQO 

J 

3 

U  1 

P  3 

0  .  '*22 1 A  5-1 7 

PA 

s 

0  .ojgagqoo 

.) 

3 

48 
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3 
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1 

3 

49 
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0.  '•223672.I 

P  A 

3 

0 .9990069-7 

3f 


DATA  OUTPUT  FOP  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 
FOR  PLATE  WITH  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  COMOITJON  IS  0.0100 
ENTIRE  SPECIMEN  MODEL  INPUT 
AC  s  4.90010000 
BETA  =  0.10000000 

ALPHA  a  20.00000000 
PSAR  *  0. 08702500 

NE  *  352 

ALPHAS  *  22.00000000 

BETAA  a  0.00289262 


J 

s 

t 

P 

s 

0.00871 700 

PA 

S 

1.00000000 

J 

s 

2 

P 

2 

0. 009(10  757 

R  A 

2 

1 .00000800 

J 

3 

3 

P 

2 

0.0092081 3 

RA 

2 

1 .00000000 

J 

3 

'l 

P 

2 

0.00058870 

PA 

2 

1.00000000 

J 

3 

5 

P 

2 

0.00047027 

R  A 

2 

1.00000000 

J 

3 

b 

P 

2 

0.01016083 

PA 

2 

1 .00000000 

J 

2 

7 

P 

2 

0.01046040 

PA 

2 

1 .00000000 

J 

2 

8 

P 

2 

0.010  75007 

R  A 

2 

1 .00000000 

J 

2 

0 

P 

2 

0.01104153 

RA 

2 

1 .ooooooon 

J 

2 

10 

P 

2 

0.01 1 S3?1 o 

PA 

2 

1 .ooooooon 

J 

2 
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P 

2 

0.01 162267 
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2 

1 .00000000 

1 

2 

12 

P 

2 

0.01 101 323 

PA 

2 

1 .onooonoo 

J 

2 

1  3 

o 

2 

0.01220380 

RA 

2 

0.09999906 

J 

2 

l'i 
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0.01 2«o«37 

PA 

2 

0. 99999990 
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2 

IS 
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2 
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2 
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p 

2 
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2 
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p 

2 
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RA 

2 
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3 

<8 

p 

2 

0.01 365663 

RA 

2 
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.1 

2 

1  o 

p 

2 
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RA 

2 

0.93999948 

J 

3 

20 

p 

2 

0.91423777 

RA 

2 

0.99999981 

J 

2 

21 

p 

2 

0.01452833 

PA 

2 

0.99999971 

I 

2 

22 

D 

2 

0.01481  89.1 

PA 

2 

0.93999954 

J 

2 

?  3 

P 

3 

0 . 0  1 5  1  0  a  u  7 

PA 

2 

0.93999930 

J 

3 

3<1 

P 

2 

0.015  10003 

PA 

2 

u.  93994994 

J 

5 

25 

O 

2 

0.01  560491. 

PA 

2 

0 .939°9«40 

J 

5 

?6 

0 

2 

0.015381 1 7 

RA 

2 

0 • 93  399  7bn 

J 

3 

?  7 

p 

2 

O.Oi 627173 

RA 

2 

0.99999643 

J 

3 

38 

p 

2 

0.01656230 

PA 

2 

3 .93999474 

J 

2 

jv 

p 

2 

0 .0 1 645?“ 7 

PA 

2 

0.93939228 

J 

3 

3<) 

p 

2 

0.11714  50  3 

RA 

2 

0.99998876 

J 

3 

31 

p 

2 

0.0  1  7434  9  <> 

RA 

2 

0.99998373 

J 

2 

32 

p 

2 

0.0  1  772457 

RA 

2 

0.93997860 

J 

2 

33 

p 

2 

0.0180151 3 

RA 

2 

0.99996653 

J 

3 

30 

p 

2 

0.018  505/0 

RA 

2 

0.93995241 

J 

2 

35 

p 

2 

0.01450,27 

RA 

2 

0.93993271 

J 

3 

3b 

p 

2 

0.018*8683 

RA 

2 

0.93990556 

J 

3 

37 

p 

2 

0.0191  77-40 

P  A 

2 

0.93086759 

.1 

3 

3* 

p 

2 

0.  J\94o7J7 

RA 

2 

0.99981567 

J 

3 

5° 

p 

2 

9 . 0 1  9  7  5  5  3 

PA 

2 

0.99974465 

J 

2 

Un 

u 

2 

0.0200  49M 

PA 

2 

0.99064795 

J 

2 

Ul 

p 

2 

0. 0?.)  5308  ’ 

PA 

2 

0. 93951644 

J 

3 

U2 

p 

2 

0..)?063n25 

PA 

2 

0.93934001 

J 

3 

<13 

0 

2 

0.0?032'®0 

PA 

2 

0.99910233 

J 

2 

U» 

p 

2 

0 . 0 2  1  2 1  l  37 

PA 

2 

0.93878427 

J 

2 

(15 

p 

2 

0.02 150133 

PA 

2 

0. 938360  39 

J 

2 

<36 

p 

2 

0.021 7025' 

PA 

2 

0.93779774 

) 

2 

<17 

D 

3 

0.  )2?r»3«7 

K  A 

2 

0.93795383 

1 

3 

48 

O 

3 

0.  122  573, < 

R  A 

2 

0.996-j/  <04 

1 

3 

40 

o 

3 

0  .  1  ?  2  3  6  <  2 

-  A 

: 

9 . 9  3  |  7  H  4  7  ,7 

/ 

3 

30 

ff 

3 

1.  -22’5  <77 

i* : 

2 

0.03514441 

data  output  for  LT»<eARize  size  errscr  uniaxial  tensile  stress 

FOR  PLATE  WITH  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  O.OOtO 
HALF  SPECIMEN  HOOEL  INPUT 


KC  s 

4. 

ooiooooo 

•ETA  * 

0. 

1QOOQOOO 

ALPHA  = 

to. 

oooooooo 

PBAR  = 

0. 

00871700 

NE  - 

A82 

ALPHA2  s 

12. 

OOOOOOOO 

BETAA  a 

0. 

03050189 

J  = 

I 

P  = 

0.0007  1  7011 

RA  s 

0.40999971 

J  a 

2 

O  - 

a. 01)000757 

RA  s 

0.99999957 

J  a 

3 

P  = 

•J.0OOJQ91  ? 

RA  s 

0.99999939 

J  S 

a 

P  2 

0 •  409  5 A  0  7  0 

RA  : 

0. 93444010 

J  a 

5 

P  2 

0.00907927 

RA  = 

0.94444871 

J  a 

6 

P  = 

O.OUU690S 

RA  s 

0.9444481 7 

J  a 

7 

P  2 

0.01  04604(1 

RA  s 

0.44444744 

J  s 

0 

P  2 

0.01't  75097 

®A  * 

0. 94444604 

J  S 

<J 

P  2 

0.01 104153 

RA  s 

0 .44444509 

J  a 

10 

□  - 

0.011 33?! 0 

RA  s 

0.99499330 

J  S 

1 1 

P  2 

0.0 H  62 ?b 7 

RA  s 

0.99499042 

.1  a 

12 

P  2 

0.01101 323 

RA  s 

0.94440779 

.1  s 

1  3 

□  - 

0.31220390 

RA  s 

0.99990369 

J  a 

1« 

P  2 

0 .0 1  ?'io<i  J7 

RA  s 

0.99997037 

J  a 

IS 

P  - 

0.01279403 

RA  s 

0.99997150 

1  a 

in 

P  2 

0.01 307550 

PA  3 

0.94946260 

J  a 

1  7 

P  2 

0. <il  336607 

RA  3 

0.44445142 

1  a 

1“ 

P  2 

0.01  3  656o  3 

RA  3 

0.99443712 

f  a 

1'4 

P  2 

0. 01304720 

W  A  s 

0.40441404 

J  a 

20 

P  2 

<>.01423777 

PA  3 

0.44464632 

.!  a 

21 

P  S 

0.0  I  <152433 

RA  = 

0.49406707 

J  a 

22 

P  2 

O.OluoiSO" 

R  A  s 

0.444«3?44 

J  a 

23 

P  = 

0.01510047 

RA  s 

0.44470647 

J  a 

?<t 

P  2 

0.015  1 0003 

RA  s 

0.43073415 

J  a 

25 

P  3 

0.0156006" 

RA  3 

0. 930bb732 

J  a 

26 

P  3 

0.0!53BI 1 7 

PA  3 

0.99450539 

J  a 

2  7 

P  3 

0.016271 73 

RA  3 

O.99O'I0534 

J  = 

2” 

P  3 

3.01656230 

PA  3 

0.99436359 

J  a 

2'' 

P  2 

0.0l6“5297 

RA  3 

0.99921595 

J  a 

io 

P  3 

<1.0  1  71  4  54  3 

RA  s 

0.99903752 

J  a 

ii 

P  3 

0.01 743400 

PA  3 

0.93092257 

J  a 

a 

P  2 

0.01772457 

RA  s 

0.99056443 

J  a 

a 

P  3 

0.01601513 

RA  s 

0.93025539 

J  a 

3*1 

P  3 

0. 01630570 

PA  3 

0.99780649 

J  a 

35 

P  2 

0.0195*4627 

RA  s 

0.94744740 

J  a 

in 

P  2 

0.0  1  9  39  6-*  3 

RA  s 

0.99692624 

J  a 

37 

P  3 

0.3  141  7 7  1" 

RA  3 

0.94630956 

J  a 

38 

P  2 

0.01*4  lb  707 

RA  3 

0.995501 10 

J  a 

55 

P  2 

0.01075953 

RA  s 

0.99472362 

J  a 

41) 

°  2 

0.02004  3  1  0 

RA  s 

0.99371655 

J  a 

4  1 

P  2 

4.020  3346  7 

?  A  3 

0.99253675 

J  a 

02 

O  2 

>.02063323 

RA  = 

0.931 15600 

J  a 

It  3 

P  2 

0 . 02*>  42  >60 

RA  3 

0.98455067 

J  a 

Alt 

P  2 

0 . 02 1 2  U  3 7 

RA  3 

O.9076«l 35 

J  a 

0.5 

P  2 

0.02  1  501*43 

RA  3 

0.98551251 

J  a 

Ob 

P  2 

0.03  1  74  25*4 

RA  3 

0.40300210 

J  a 

a  7 

P  3 

0 . 0 2 2 ■> -  30  7 

RA  3 

0.9901  031  0 

J  a 

it* 

P  3 

0.422  <730s 

RA  3 

0.97076349 

J  a 

49 

P  3 

0  .  <*  32  6r»  130 

V  A  3 

0.9729251* 

OATA  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 
FOR  PLATE  WITH  NOTCH 

INHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  0.0040 
HALF  SPECIMEN  MODEL  INPUT 
KC  s  0.00100000 

SETA  *  0.10000000 

ALPHA  =  10.00000000 

PBAR  :  0 . 03485000 


NE 

S 

482 

ALPHAS 

s 

12. 

00000000 

SETAA 

3 

0. 

03074947 

J 

3 

1 

P  3 

0.0)871 70  0 

64 

3 

0.99999973 

J 

3 

2 

P  3 

0.0O9-J0  7S7 

68 

s 

0.999999b0 

J 

3 

5 
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data  OUTPUT  FOR  LINEARfZE  SIZE  EFFECT  UN t AXIAL  TENSILE  STRESS 
FOR  PLATE  WITH  NOTCH 

WHERE  THE  DISPLACEMENT  0OUNOARY  CONOITJON  13  O.OIOO 
WALF  SPECIMEN  M006L  IWPUT 
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OAT*  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAKUL  TENSILE  STRESS 
FOR  PLATE  Y*IT«  NOTCH 

WHERE  THE  DISPLACEMENT  boundary  conottion  r$  0.01 00 
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DATA  OUTPUT  FOR  tlN£A«I2£  SIZE  EFFECT  UNIAXIAL  TENSILE  STRESS 
FOR  PLATE  WITHOUT  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  O.OOIO 
ENTtRE  SPECIMEN  MODEL  INPUT 
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Data  OUTPUT  FOR  LINEARIZE  SIZE  EFFECT  UNIAXIAL  tensile  stress 
FOR  PLATE  WITHOUT  NOTCH 

WHERE  THE  DISPLACEMENT  BOUNDARY  CONDITION  IS  O.OIOO 
ENTIRE  SPECIMEN  MODEL  INPUT 
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BETA  x  0.10000000 
ALPHA  x  10.00000000 
PBAR  x  0.00491000 
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COMPUTER  PROGRAMS  FOR  CONVERTING 

FINITE  ELEMENT  (NIKE2D)  PROGRAM  OUTPUT  TO  BE 

STATISTICAL  RELIABILITY  COMPUTER  PROGRAM  INPUT 

. . . . . . . 

c 

c  PROGRAM:  CONVERT  MTKE3D  OUTPUT  ro  BE  3TAT1S1CAL  INPUT 

C 

C  FOR  1FT.  INTEGRATION  ENTIRE  SPECIMEN  MODEL  ONLY 


■  e 


XW<M>=.  02*4375 
lt><>  CONTINUE 

00  no  1**317.328 
XL  (N)  * .1658334 
X#(m>=.02457S 
liu  CONTINUE 

DO  III  N*329.360 
XL  (M>  =.0995 
XMfM)*.02487S 

1 1 1  CONTINUE 

DO  422  N=J6i.42a 
XL(M) =.024875 
XW(N) *.024875 

112  CONTINUE 

00  113  *=425,456 
XL tN) =.0905 
XW(M/ =.024075 

113  CONTINUE 

DO  114  H*457»  488 
XL (M)=. 1658334 
XN(M) =.024875 
It 'I  CONTINUE 

00  1  IS  4=464.576 
XL <»> 3.0*4875 
xmN)  s.024875 

115  CONTINUE 

DO  1 16  *=577,592 
XL(M)=. 024875 
XN<  Mi  =  .024875 

1 1 6  CONTINUE 

DO  10  Isl.NP 
00  13  7=1.10 

Ms  19  »  ( I- 1 )  ♦  j 
IFlH.Nt.NEi  TM(N 

IF(J.£*.10)  then 

*£40(5.45)  SICY.SIGZ.SIlX.SIGYZ 

ELSE 

RCA0(5.35)  SIGy.SIGZ.SIGX.SIGyz 
ENO  IF 

ELSE 

BEAD15.30)  SIGV.S1GZ.SIGX.SIGVZ 
CO  TO  55 
ENO  IF 

WRIT  EC  6, 40  )  N.Xt-(M)  ,**<*) .SIGY.SIGZ.SIGX.SIGYZ 
IS  CONTINUE 
10  CONTINUE 

55  **176(4.40)  N, XL(N) ,X*(M1 .STCY.SIGZ.SIGy.SIGYZ 
30  FORMAT ( iSX.ftEll .0) 

35  FORHAT( 18X, 4E 1 1.0.////) 

40  FORMAT (.3%,  I5.2F1  I  .7.4E12.4) 

»5  FORMAT ( 1 3 < , 4t1 1.0, ////////////) 

50  FORMAT  (  IX.  ’ELEMENT  N0.’,4X,  'LMOl,  'W.7X,  ’SIGY’.a*.  ’SIGZ 
•  ,8X, 'SIGX',7X. ' SIGVZ ' » Z) 

STOP 


o  o  o  o  o  o  o 


PROSRAH:  CONVERT  NIKE20  OUTPUT  TO  »£  STaTISI 
FOR  4 pr.  INTEGRATION  half  specimen 


c 

c 

. . . . . . . 

c 

Dimension  %L(»925l.*Wtt*2S) 
syte  a<30> 

TrPfi., 'ENTER  FILE  CONTAINING  NIKE?!)  OUTPUT' 
REAO(5,5«0)  A 
Son  F0RNAT(3iUi  i 

OPEN  ( UN  l  T  aS »  TUCA,  STAIUSs'OLP' ) 

OPEN(lJNI  T s6.  FILE5’  M.  0&  I  '  »  S  t  A  TUSs  •  NEW  ) 

MR  1 TE  ( b,  50  > 

NFrfftP 

AE«=NE»4 
NPsNE4/40»t 
00  tOQ  Nsi , 100 
XL ( M ) 3 .0*975 
XW(M)=. 02850 

1 00  CONTINUE 
DO  101  M=101 .220 

XL  ( f#> s .024875 
X«NMl  =  .028500 

101  CONTINUE 

00  102  *=211,540 
xuN'  =  .  01 24375 
XW<M  >  = . 0285  000 

102  CONTINUE 
00  1^3  1*3341,460 

XL (N) 3.024875 
XNfNls.028500 

103  CONTINUE 
00  1  04.  M 3661.760 

xL(M)s.0497S 
xW(M) s.02850 
l  oa  CONTINUE 

00  105  Ms  76 1 , 920 
XV fM I s. 04975 
XN(M> 3.0124375 
105  CONTINUE 

00  lOfe  *=921 ,1112 
Xt(Mls.02a875 
XW(MI  =  .0 12*375 
I  Oo  CONTINUE 

00  107  M=1 t 13, 1568 
XL  ( Ml  s ,  0 1 243  75 
XW(M) s.OI 24375 
107  CONTINUE 

00  1*® *=1364, 1560 
XLTMTs. 02*875 
X*(  n  )  s  .0124375 

HIP  continue 

00  109  Ms  150, 1720 


XLlM) =.04975 
XW(M)=.012437S 
10*  CONTINUE 

00  110  M=1 721, 1896 
XUm  =  .012«75 
XW(M)=.0l2437S 

110  CONTINUE 

00  111  *=1897.1928 
XL (M) =.01 24375 
XW(M>=. 0124375 

1 1 1  CONTINUE 

00  10  1=1, NP 
00  15  3=1,10 
00  20  <=1,4- 

Ma40*(l-I  J+4»(J-n*K 
IF  (M.NE.NEt)  THEN 
I F  <  |C.  EQ  .«  )  THEN 

fFtj.Eo.io)  then 

READ (5, 45)  SIGY,3IG7,SIGX,SIGYZ 

else 

READ (5,35)  SIGT.5IGZ, SICX, SIGYZ 
END  IF 
ELSE 

READ (5, 30 )  Sl&Y.SIGZ.SIGx.siGYZ 
ENO  IF 
ELSE 

READ fS ,30 )  SIGY.SIGZ.STGX.STGYZ 
GO  TO  55 
ENO  IF 

*R|TE(6.40j  *,XL(M) ,XWCN) ,SIST,STGZ,SIGX,31GYZ 
20  CONTINUE 

IS  CONTINUE 
10  CONTINUE 

55  V»fiIT£<6,40)  M,XHM).XWfN),5lG»,SIGZ,SrGX,SI5YZ 
30  F0«MAT(jaX,4E11  .0) 

35  FORMAT ( tan, 4E 11 .0. O 
*0  FORMAT (3X. 15. 2F11  .7.4E12.4) 

«5  FORNATf 18X,4E11 .0, /////////) 

50  FORMATf tx, ’ELEMENT  NO. ' . 4X, 'L ' , I  OX, >W> ,7X, ’ SIGT • ,ax , • SI 
*,ax,lSICx•,7x,,SIGYZ•,/) 

STOP 

ENO 


IK6P 


non 


PROGRAM:  CONVERT  NIKE2D  OUTPUT  TO  BE  STATISTICAL  I*»PUT 

TOR  iPT.  INTEGRATION  MALE  SPECIMEN  MODEL  ONLY 
C 

. . . . . 

c 

c»»**»*****«******.«****»» *»*»»•**••*** »••»»•*•* 

c 

DIMENSION  XL ( 462) » XW( 492) 

BYTE  A(iO> 

rrPE»,'FIL£  CONTAJNINE  NIKE 2D  OAT A  OUTPUT* 

READ(5,SOO )  A 
S00  FORMAT (JOAl) 

OPENfONITsS,  FILE*A.  STATUS=*OLO* ) 

OPENlUNlTsb,  FILE* 1 aB.DAT' »  STATU6* *NEM* ) 

WRITE (6, SO) 

NE=4« 

NPsNC/1 0*1 
00  140  Ms! #25 
XLfM)=.099S 
Xy»(M)  =  .0570 

100  CONTINUE 
DO  101  M=?i,SS 

XUM)s.  04973 
XW(M)=.0S7 

10 1  CONTINUE 
DO  102  M=56. 13S 

XL fM)=. 02*875 
Xw(Mi=.4570 

1 02  CONTINUE 
00  143  M=136. IDS 

XUM>  =  . 04475 
XW(N ) = .0570 

1 03  CONTINUE 
DO  104  M= I  66 # 190 

XL  t  M )  =  . 0995 
XW(M) =.0570 
10U  CONTINUE 

DO  lOS  Ms  191, 230 
XL (M) =.0995 
XW(M) =.024875 

105  CONTINUE 
00  lOo  M=23l .278 

XL(N)s. 04973 
XW(M)s. 024D75 

106  CONTINUE 
00  107  Ms279,l42 

<L(M)=. 02*675 
XW<  M ) s ,024075 

107  CONTINUE 
00  1 OF  4=343. 140 

XLfN)=. 04975 
XWCM) =.024875 
106  CONTINUE 

»0  100  Ms 341 . 43 0 
XL'M) = .0495 


80 


XX  CM  >  =  .024*75' 

10*  CONTINUE 

00  110  *=431,474 
Xl(M) S.02487S 
XW(N)=. 024*75  . 

110  CONTINUE 

DO  111  M-475. 482 
XU (M) =.024*75 
XWfMIs. 024*75 

111  CONTINUE 

00  10  1*1, NP 
00  15  J=1.10 

N=1 *J 
IF(N.NE.NE)  THEN 

IFU. £0.101  THEN 

REA0<5,«S)  SIGY,SIGZ,SIGX,SIGYZ 

ELSE 

R£AD(5,3S)  SIGY ,SICZ,SICx,SIGY2 
CNO  IF 

ELSE 

R£AD(5 ,301  SIGY,SIGZ,SIGX,SIGYZ 
GO  TO  55 
CNO  IF 

WRITE! 6, 40)  M,XL(M),XM(M>,SIGY,SIGZ,3IGX,31GY; 

15  CONTINUE 
10  CONTINUE 

55  XRITEC6.40)  M,XLCN),XW(H1,3ICY,SISZ,SIGX,SIGYZ. 

3o  FoRM4r(i*x.4«n.o) 

35  F08M*T(18X,461I  .0,////) 

00  F0RMAT(3X. 15, 2F11. 7,4612. 4) 

05  FORMAT! iaX,4EH  .0, ////////////) 

50  FORNATfIX. 'ELEMENT  NO. • ,4X, 'L' , I  OX, , 7x , 1 SIGY • ,*X , • SIGZ 
•,8X. 'SIGX' ,7X, 'SIGYZ' ,/) 
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